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and adhesion for rough surfaces, where the power spectrum enters as an important 
input. 
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1. Introduction 



Surface roughness has an enormous influence on many important physical phenomena 
such as contact mechanics, sealing, adhesion and friction. Thus, for example, 
experiments have shown that already a substrate with a root-mean-square (rms) 
roughness of order ~ 1 /im can completely remove the adhesion between a rubber 
ball and a substrate, while nanoscale roughness will remove the adhesion between most 
hard solids, e.g., metals and minerals; in other words, roughness is the main reason why 
adhesion is usually not observed in most macroscopic phenomena. Similarly, rubber 
friction on most surfaces of practical interest, e.g., road surfaces, is mainly due to the 
roughness-induced pulsating forces which act on the rubber surface as it slides over the 
substrate asperities. 




Figure 1. (a) Micrometer sized cantilever beam, (b) If the beam is too long or too thin 
the minimum free energy state corresponds to the beam partly bound to the substrate. 
Surface roughness lowers the binding energy (per unit area) and hence stabilizes the 
non-bonded state in (a). 

Let us illustrate the importance of surface roughness with three modern 
applications. At present there is a strong effort to produce small mechanical devices, e.g., 
micromotors. The largest problem in the development of such devices is the adhesion 
and, during sliding, the friction and wear between the contacting surfaces [1]. As an 
example, in Fig. 1 we show the simplest possible micro device, namely a micrometer 
cantilever beam. (Suspended micromachined structures such as plates and beams are 
commonly used in the manufacturing of pressure and accelerator sensors.) If the beam 
is sufficiently long or thin the free beam state in (a) will be unstable, and the bound 
state in (b) will correspond to the minimum free energy state [2]. Roughly speaking, the 
state (b) is stable if the binding energy to the substrate is higher than the elastic energy 
stored in the bent beam. The binding energy to the substrate can be strongly reduced 
by introducing (or increasing) the surface roughness on the substrate (see Sec. 7.1). In 
addition, if the surfaces are covered by appropriate monolayer films, these surfaces can 
be made hydrophobic thus eliminating the possibility of formation of (water) capillary 
bridges. 

A second application is the formation of hydrophobic coatings on surfaces by 
creating the appropriate type of surface roughness [3] . This amounts to copying Nature 
where many plant surfaces are found to be highly hydrophobic (Fig. 2) as a result of the 
formation of special types of surface roughness (Fig. 3). The surface roughness allows air 
to be trapped between the liquid and the substrate, while the liquid is suspended on the 
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Figure 2. A water droplet on a superhydrophobic surface: The droplet touches the 
leaf only in a few points and forms a ball. It completely rolls off at the slightest 
declination [3]. Reproduced with permission. 




Figure 3. A leaf surface with roughness on several length scales optimized (via 
natural selection) for hydrophobicity and self-cleaning. Through the combination of 
microstructure (cells) and nanostructure (wax crystals) the water contact angle 9q is 
maximized [3] . Reproduced with permission. 

tips of the asperities. Since the area of real liquid-substrate contact is highly reduced, 
the contact angle of the drop is determined almost solely by the surface tension of the 
liquid, leading to a very large contact angle. New commercial products based on this 
"Lotus effect" [4], such as self-cleaning paints and glass windows, have been produced. 

Finally, we mention the effect of surface roughness on direct wafer bonding [5]. 
Wafer bonding at room temperature is due to relatively weak interatomic attraction 
forces, e.g., the van der Waals interaction or hydrogen bonding, giving (for perfectly flat 
surfaces) an interfacial binding energy of order 6 meV / A . The wafer surface roughness is 
the most critical parameter determining the strength of the wafer bonding. In particular, 
when the surface roughness exceeds a critical value, the wafers will not bind at all, in 
agreement with the theory presented in Sec. 7.1. Primary grade polished silicon wafer 
surfaces have rms roughness of order ~ 0.1 nm when measured over a 10 x 10 /im surface 
area, and such surfaces bind spontaneously. However, when the surface roughness 
amplitude is of order 1 nm the surfaces either bind (slowly) when squeezed together 
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at high enough pressure, or they do not bind at all depending on the detailed nature of 
the surface roughness power spectra. 

Surfaces with "ideal" roughness, e.g., prepared by fracture or by some growth 
process, have been studied intensively for many years [6, 7, 8, 9]. However, much 
less information has been available for more common surfaces of engineering interest. 
In this article we discuss the nature of the power spectra of some surfaces of practical 
importance. As illustrations we discuss contact mechanics, sealing, rubber friction and 
adhesion. 

This paper is organized as follows: In Sec. 2 we define the surface roughness power 
spectrum C(q), and discuss some of its properties. In Sec. 3 we present power spectra 
deduced from surface topography measurements for different surfaces of technological 
importance. Sec. 4 illustrates how the surface roughness power spectrum determines 
the contact area between elastic and viscoelastic solids. In Sec. 5 we briefly consider the 
influence of surface roughness on sealing. In Sec. 6 we discuss rubber friction on rough 
substrates, and show some results for how the friction depends on the surface roughness 
power spectrum. We also discuss the influence of tire-road polishing and discuss rubber 
friction on wet road surfaces. In Sec. 7 we discuss adhesion between rough surfaces 
and present two applications which illustrate the importance of surface roughness for 
adhesion in technology and biology. Sec. 8 contains our summary and an outlook. 

2. Surface roughness power spectrum: definition and general properties 

The influence of roughness on the adhesional and frictional properties introduced above 
is mainly determined by the surface roughness power spectrum C(q) (or power spectral 
density) defined by [10] 



Here x = (x, y) and z = h(x) is the substrate height measured from the average surface 
plane, defined so that (h) = 0. The (. . .) stands for ensemble averaging, i.e., averaging 
over a collection of different surfaces with identical statistical properties. We have 
assumed that the statistical properties of the substrate are translationally invariant, so 
that the correlation (h(x + x )h(x )) does not depend on the choice of x Q , but only 
on the in-plane distance vector x. Alternatively the height correlation can be defined 
through a spatial average over x : 



The latter definition is quite popular in the engineers' community [10], and it is closer to 
the experimental procedure to measure the power spectrum. We notice, however, that 
the lateral size L of the surface area is always finite in any experimental measurement 
or numerical calculation. In some special the ones discussed in Sec. 3.1, the 

spatial average is not equivalent to the ensemble average for any practical choice of L. 
Finally we point out that the power spectrum can also be defined through the square 




(1) 
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modulus of the Fourier transform of h; this is indeed the conventional approach adopted 
in the context of signal theory and electronics, and it explains the origin of the name 
"power spectral density". Nonetheless such approach is equivalent to Eq. (1) because 
of the Wiener-Khintchine theorem (see Appendix C and particularly Eq. (C.6)). 

Together with translational invariance, in (1) we have also assumed that the 
statistical properties of the substrate are isotropic, so that C(q) only depend on the 
magnitude q = \q\ of the wave vector q. Note that from (1) follows 

(h(x)h(o)) = J d 2 q C{q)e iqx 

so that the root-mean-square roughness amplitude a = (h 2 ) 1 / 2 is determined by 

d 2 qC(q) = 2nj o dq qC{q) (2) 

In reality, there will always be an upper limit and a lower limit to the g-integral in 
(2). Thus, the largest possible wave vector will be of order 27r/a, where a is a short 
wavelength cutoff corresponding perhaps to some lattice constant, while the smallest 
possible wave vector is of order 2n/L where L is the linear size of the surface. In 
general, one may define a root-mean-square roughness amplitude which depends on the 
range of roughness (?o,?i) involved in the integral in (2): 

(h 2 }(q , qi ) = 2n f 9l dqqC(q) (3) 
J go 

For a randomly rough surface, when h(x) is a Gaussian random variable, the 
statistical properties of the surface are completely defined by the power spectrum C(q) 
(see Appendix A). In this case the height probability distribution 

P h = {S[h-h(x)]) 

will be a Gaussian 

p _ 1 c -hW 
h (2vr)VV 

The height distribution of many natural surfaces, e.g., surfaces prepared by fracture, 
or surfaces prepared by bombardment by small particles (e.g., sand blasting or ion 
sputtering) are usually nearly Gaussian. On the other hand, rough surfaces, e.g., 
prepared by fracture, and then (slightly) polished, display a non-symmetric height 
distribution (i.e., no symmetry as h — > —h) since the asperity tops have been polished 
more than the bottom of the valleys. These surfaces, typically of considerable practical 
importance-see below, have non-Gaussian height distribution. For such surfaces it is 
interesting to study the top, Ct, and the bottom, C B , power spectra defined by 

C T (q) = J^J^ (h T (x)h T (o))e-^ x (4a) 

C B (q) = J^yJ d2x (h B {x)h B {o))e-^ x (46) 

where Ht{x) = h(x) for h > and zero otherwise, while h B (x) = h(x) for h < and 
zero otherwise, the latter being 'rectified' profiles, see Fig. 4. It is easy to show that 
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Figure 4. The surface profile h(x) is decomposed into a top hx{x) and a bottom 
}ib{x) profile. 



C ~ Ct + Cb- It is also clear by symmetry that for a surface prepared by fracture, 
Cr(q) = Csiq), since what is top on one of the cracked block surfaces is the bottom on 
the other (opposite) crack surface, and vice versa, see Fig. 5. However, if the cracked 
surface is slightly polished then, because the contact pressure with the polishing object 
(e.g., sand paper) is highest at the asperity tops, the polishing at the asperity tops will 
be stronger than at the valley bottoms, and thus Cb > Ct- If and % are the fraction 
of the nominal surface area (i.e., the surface area projected on the xy-plane) where h > 
and h < 0, respectively, with ut + % = 1, then we may also define C T {q) = Ct/ut 
and C* B = C B /n B - In general, 0.5 and for surfaces prepared by fracture 

ut = ns = 0.5. Roughly speaking, C T would be the power spectrum resulting if the 
actual bottom profile (for h < 0) was replaced by a mirrored top profile (for h > 0). A 
similar statement holds for C B . 

Many surfaces tend to be nearly self-afline fractal. A self-afline fractal surface has 
the property that if part of the surface is magnified, with a magnification which in general 
is appropriately different in the perpendicular direction to the surface as compared to 
the lateral directions, then the surface "looks the same", i.e., the statistical properties 
of the surface are invariant under the scale transformation (see Appendix B). For a 
self-afline surface the power spectrum has the power-law behaviour 



where the Hurst exponent H is related to the fractal dimension Df of the surface via 
[6] H — 3 — Df. Of course, for real surfaces this relation only holds in some finite wave 
vector region q < q < qi, and in a typical case C(q) has the form shown in Fig. 6. Note 
that in many cases there is a roll-off wavelength q below which C(q) is approximately 
constant. We will discuss this point further below. 

Finally, note that while the root-mean-square roughness usually is dominated by the 
longest wavelength surface roughness components, higher order moments of the power 
spectra such as the average slope or the average surface curvature are dominated by 
the shorter wavelength components. For example, assuming a self afline fractal surface, 
Eq. (3) gives 



if qi/qo ^> 1. However, the average slope and the average curvature have additional 
factors of q 2 and q 4 , respectively, in the integrand of the ^-integral, and these quantities 



C{q) ~ q 



2(H+1) 
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cracked 
block 





Figure 5. Rough surfaces prepared by crack propagation have surface roughness with 
statistical properties which must be invariant under the replacement of h — > — h. This 
follows from the fact that what is a valley on one of the crack surfaces (say the left) is 
an asperity with respect to the other crack surface (right). Thus the top and bottom 
power spectra must obey Ct{<i) = Cb(q)- 



are therefore dominated by the large q (i.e., short wavelength) surface roughness 
components (see Appendix C). 

3. Surface roughness power spectrum: experimental results 

In this section we present power spectra for different surfaces of practical importance. 
The power spectra have been calculated using Eqs. (1), (4a) and (4b) (see Appendix 
D), where the height profile h(x) has been measured using either optical methods^ or 
by Atomic Force Microscopy [9]. 

3.1. Surfaces produced by crack propagation 

Fig. 7 shows the power spectra C(q) for three freshly cleaved stone surfaces, namely a 
granite and two basalt stone surfaces. Here, and in what follows, we show the power 
spectra on a log-log scale, where log stands for the base 10 logarithm. Note that the 
granite and basalt surfaces, in spite of the rather different mineral microstructure (see 
below), give identical power spectra within the accuracy of the measurement. It has been 

| Many methods for measurement of surface topography have been developed see, e.g., 
http://www.michmet.com, http://www.solarius-inc.com or http://www.schmitt-ind.com for experi- 
mental equipment and some illustrative results. 
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Figure 6. Surface roughness power spectrum of a surface which is a self-affine fractal 
for 50 < q < q%. The long-distance roll-off wave vector q and the short distance 
cut-off wave vector q\ depend on the system under consideration. The slope of the 
logC — log q plot for q > q determines the fractal exponent of the surface (see text). 
The lateral size L of the available surface region determines the smallest possible wave 
vector qL = 2tt/L. 




Figure 7. The surface roughness power spectra for two freshly cleaved basalt surfaces 
and a fresh granite surface. 
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stated (see, e.g., Ref. [11]) that surfaces produced by crack propagation have self affine 
fractal structure with the universal fractal dimension Df ps 2.2. However, our measured 
logC — logg relations are not perfectly straight lines, i.e., the surfaces in the length- 
scale range studied cannot be accurately described as self-affine fractals. Moreover the 
average slope of the curves in Fig. 7 correspond to the fractal dimension D{ ps 2 rather 
than 2.2. 

The similarity of the power spectra for the basalt and granite surfaces in Fig. 7 is 
striking. Granite and basalt both result from magma and have a similar composition, 
consisting mainly of minerals from the silicate group. However, granite originates from 
magma which was trapped deep in the crust, taking a very long time to cool down and 
crystallize into solid rock. As a result granite is a coarse-textured rock where individual 
mineral grains are easily visible. Basalt, on the other hand, results from fast cooling of 
magma from, e.g., volcanic eruptions, is fine grained, and nearly impossible to resolve 
into the individual mineral grains without magnification. In spite of these differences, 
the surface roughness power spectra of freshly cleaved surfaces are nearly identical. This 
may indicate some kind of universal power spectrum for surfaces resulting from cleaving 
of mineral stones of different types, a point that could deserve further investigation. 

Note that there is no roll-off region for these fracture-produced surfaces, whose 
behaviour appears fractal-like up to the longest length scale studied, determined by the 
lateral size L of the surfaces (or of the regions experimentally studied, of the order of 
1 cm in our case), i.e., with reference to Fig. 6, go — Ql = 27r/L. One consequence of 
this is that the rms-roughness amplitude is determined mainly by the A ~ L wavelength 
fluctuations of the surface height, and will therefore depend on the size L of the surface. 
Furthermore, the height distribution Ph obtained for any given realization of the rough 
surface will not be Gaussian, but will exhibit random fluctuations as compared to other 
realizations (see Fig. 8, which illustrate this point for the three stone surfaces discussed 
above). However, the ensemble averaged height distribution (not shown) should be 
Gaussian or nearly Gaussian. Thus, when there is no roll-off region in the measured 
power spectra, averaging over the surface area is not identical to ensemble averaging. 
However, when there is a roll-off wave vector q = 2ir/X , and if the surface is studied 
over a region with the lateral size L ^> Ao, ensemble averaging and averaging over the 
surface area L x L will give identical results for Ph, and the rms-roughness amplitude 
will be independent of L for L ^> A . 

3.2. Polished crack surfaces 

In the past, cobble stones, made of granite or basalt, were frequently used for road 
surface pavements. However, these surfaces do not exhibit good frictional properties 
against rubber. In particular, with increasing time of use, the cobble stone surfaces 
become polished by slipping tires (see Sec. 6.2), and that polishing results in a reduced 



On the nature of surface roughness with applications 



11 
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Figure 8. The height distribution Ph for two freshly cleaved (cobble stone) basalt 
surfaces and a fresh granite surface. Note the random non-Gaussian nature of the 
height profiles. 




Figure 9. The surface roughness power spectra C(q) for two freshly cleaved cobble 
stone (basalt) surfaces, and for a wear-polished (used) surface. 
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Figure 10. The top and the bottom Cg surface roughness power spectra C(q) for 
a used cobble stone (basalt) surface. 

rubber-road friction, even during dry driving conditions§. Fig. 9 illustrate this polishing 
effect. It shows the power spectrum of a strongly used (basalt) cobble stone, and of two 
freshly cleaved surfaces (from Fig. 7), from the same cobble stone. At long wavelength 
the power spectrum of the strongly used surface is nearly one decade smaller than 
that of the freshly prepared surfaces. The effect of polishing is further emphasized by 
calculating the top and bottom power spectra, C£ and C* B , as shown in Fig. 10. The 
top power spectrum of the worn surface is a factor ~ 30 times smaller than the bottom 
spectrum for all wave vectors. As anticipated, the asymmetry arises from the higher 
polishing of asperities relative to valleys. It is of course crucial to take this polishing 
effect into consideration when designing road pavements, and we will discuss this point 
further in Sec. 6.2. 

3.3. Surfaces with long-distance roll-off 

As pointed out above, surfaces prepared by fracture have no natural long-distance cut- 
off and the rms roughness amplitude increases continuously and without limit as the 
probed surface area increases. This is similar to Brownian motion where the mean square 
displacement increases without limit (as ~ t 1 ^ 2 ) as the time t increases. However, most 
surfaces of engineering interest have a long distance cut-off or roll-off wavelength Ao 
corresponding to a wave vector go — 27r/Ao, as shown in Fig. 6. For example, if a 
flat surface is sand blasted for some time the resulting rough surface will have a long 

§ Cobble stones made from porphyry (a volcanic rock), thanks to the difference in hardness of the 
minerals it contains (see Sec. 6.2), have better polishing properties than cobble stones made of granite 
and basalt. See, e.g., http://www.bourgetbros.com. 
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Figure 11. The surface roughness power spectra for a fresh granite surface and a 
fresh sintered corundum surface. 

distance roll-off length, which increases with the time of sand blasting. Similarly, if 
atoms or particles are deposited on an initially flat surface the resulting rough surface 
will have a roll-off wavelength which increases with the deposition time, as has been 
studied in detail in recent growth models. Another way to produce a surface with 
a long-distance roll-off wavelength is to prepare the solid as a conglomerate of small 
particles. A nominally flat surface of such a solid has still roughness on length scales 
shorter than the diameter of the particles, which therefore may act as a long distance 
roll-off wavelength. We illustrate this here with a solid produced by sintering together 
corundum particles at high temperature and pressure (Fig. 11), and for a sandpaper 
surface (Fig. 13). For both surfaces the height distribution P h is smooth (see Figs. 12 
and 14), since averaging over a surface area with lateral size L ^> A is equivalent to 
ensemble averaging. 

The sandpaper surface in Fig. 13 was studied using the AFM at two different 
resolutions over square areas 20 x 20 /im and 100 x 100 /im as indicated by the two 
different lines in Fig. 13. The height distribution (and hence also the rms-roughness 
amplitude) calculated from these two different measurements over different surface areas, 
see Fig. 14, are nearly identical, as expected when L is larger than the roll-off length Ao- 

3.4- Road surfaces 

Asphalt and concrete road pavements have nearly perfect self-affine fractal power 
spectra, with a very well-defined roll-off wave vector q = 27r/Ao of order 1000 m -1 , 
corresponding to A ~ 1 cm, which reflects the largest stone particles used in the asphalt. 
This is illustrated in Fig. 15 for two different asphalt pavements. From the slope of the 
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Figure 12. The height distribution Ph as a function of the height h for a sintered 
corundum surface. 




log q (1/m) 

Figure 13. The surface roughness power spectra C(q) for a sandpaper surface. The 
two curves are based on the height profiles measured with an AFM at two different 
spatial resolution over 20 x 20 (red) and 100 x 100 /im (green) square areas. 
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Figure 14. The surface roughness height probability distribution Ph for a sandpaper 
surface. The two curves are based on the height profiles measured with an AFM at 
two different spatial resolution over 20 x 20 (red) and 100 x 100 /jm (green) square 
areas. 

curves for q > qo one can deduce the fractal dimension D{ w 2.2, which is typical for 
asphalt and concrete road surfaces. The height distributions of the two asphalt surfaces 
are shown in Fig. 16. Note that the rms roughness amplitude of surface 2 is nearly twice 
as high as for surface 1. Nevertheless, the tire-rubber friction is slightly higher on the 
road surface 1 because it has slightly larger power spectra for most g-values in Fig. 15. 
Thus there is in general no direct correlation between the rms-roughness amplitude and 
the rubber friction on road surfaces, as will be further discussed in Sec. 6. 

Many attempts have been made to relate rubber friction on road surfaces to the so 
called "sand filling number" . The sand filling number is the amount of very fine-grained 
sand needed to fill out all the road surface cavities in a given surface area. However, 
no correlation between the sand filling number and rubber friction on dry road surfaces 
has been found [12]. In the light of modern rubber friction theories, this result is not 
unexpected since the rubber friction depends on the power spectrum for all wave vectors, 
while only the long-wavelength components contribute appreciably to the sand filling 
number. Thus, for example, of the two asphalt surfaces in Fig. 15, the sand filling factor 
of road 2 is nearly twice as large as for road 1, but the rubber friction is slightly higher 
on road surface 1. 
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Figure 15. The surface roughness power spectra C(q) for two asphalt road surfaces. 
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Figure 16. The height distribution Ph for two different asphalt road surfaces. 



3. 5. Other surfaces of practical interest 

Finally, let us consider two other surfaces of practical importance. Fig. 17 shows 
the power spectrum of a plexiglas surface measured (using an AFM) at two different 
resolutions over two different surface areas 20 x 20 /iin and 100 x 100 fxm wide. This 
surface does not exhibit a roll-off wave vector in the studied wave vector range and 
the height distributions deduced from the two different surface areas differ strongly, 
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Figure 17. The surface roughness power spectra C(q) for a plexiglas surface. The 
two curves are based on the height profiles measured with an AFM at two different 
spatial resolution over 20 x 20 and 100 x 100 fim square areas. 
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Figure 18. The surface roughness height probability distribution P/, for a plexiglas 
surface. The two curves are based on the height profiles measured with an AFM at 
two different spatial resolution over 20 x 20 and 100 x 100 /im square areas. 

see Fig. 18, with the rms roughness amplitude being more than twice as large for the 
measurement over the larger surface area. The logC-logg relation in Fig. 17 is a nearly 
perfect straight line, and the slope correspond to the fractal dimension Df « 2.4. 
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Figure 19. The top Cj, and the bottom C* B surface roughness power spectra for a 
polymer film spin coated on a very smooth substrate and dried. From Ref. [13]. 

Fig. 19 shows the top and bottom power spectra and C* B for a polymer film 
which was spin coated and dried on a flat silicon surface. In this case there is a roll-off 
wavelength Ao ~ 1 /um which is probably related to the average polymer film thickness. 
Note that within the accuracy of the experiment, = Cg, proving qualitatively 
that the short-wavelength roughness in the large valleys is very similar to the short- 
wavelength roughness on the large asperities. The logC — logg relation in Fig. 19 for 
logg > 6.6 is a nearly perfect straight line, and the slope corresponds to a roughness 
exponent H = 1.6. This is larger than unity, which implies that the surface cannot be 
described as fractal at any scale. 

4. Contact mechanics 

Practically all macroscopic bodies have surfaces with roughness on many different length 
scales. When two bodies with nominally flat surfaces are brought in contact, real 
(atomic) contact will only occur in small randomly distributed areas, and the area of 
real contact is usually an extremely small fraction of the nominal contact area. We can 
visualize the contact regions as small areas where asperities from one solid are squeezed 
against asperities of the other solid; depending on the conditions the asperities may 
deform elastically or plastically. 

How large is the area of real contact between a solid block and a substrate? This 
fundamental question has extremely important practical implications. For example, it 
determines the contact resistivity and the heat transfer between the solids. It is also of 
direct importance for wear and sliding friction [14], e.g., the rubber friction between a 
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tire and a road surface, and has a major influence on the adhesive force between two 
solid blocks in direct contact. 

Contact mechanics has a long history. The first study was presented by Hertz [15]. 
He gave the solution for the frictionless normal contact of two elastic bodies of quadratic 
profile. He found that the area of real contact AA varies nonlinearly with the load or 
squeezing force: AA oc F$/ 3 . In 1957 Archard [16] applied the Hertz solution to the 
contact between rough surfaces and showed that for a simple fractal-like model, where 
small spherical bumps (or asperities) where distributed on top of larger spherical bumps 
and so on, the area of real contact varies nearly linearly with F N . A similar conclusion 
was reached by Greenwood and Williamson [17, 18, 19] who again assumed asperities 
with spherical summits (of identical radius) with a Gaussian distribution of heights, as 
sketched in Fig. 20(b). A more general contact mechanics theory has been developed 
by Bush et al. [20] They approximated the summits by paraboloids and applied the 
classical Hertzian solution for their deformation. The height distribution was described 
by a random process, and they found that at low squeezing force Fn the area of real 
contact increases linearly with F N . 

Fig. 21 shows the contact between two solids at increasing magnification (. At low 
magnification (( — 1) it looks as if complete contact occurs between the solids at many 
macro asperity contact regions. When the magnification is increased and smaller length 
scale roughness is detected, it can be observed that only partial contact occurs at the 
asperities. In fact, if there were no short distance cut-off the true contact area would 
eventually vanish. In reality, a short distance cut-off always exists, e.g., the interatomic 
distance. In many cases the local pressure in the contact regions at the asperities may 
become so high that the material yields plastically before reaching the atomic dimension. 
In these cases the size of the real contact area will be determined mainly by the yield 
stress of the solid. 

4-1. Elastic contact mechanics 

From contact mechanics (see, e.g., Ref. [19]) it is known that in the frictionless contact 
of elastic solids with rough surfaces the contact stresses depend only upon the shape 
of the gap between them before loading. Thus, without loss of generality, the actual 
system may then be replaced by a flat elastic surface [elastic modulus E and Poisson 
ratio u, related to the original quantities via (1 — v 2 )/E = (1 — v\)/E\ + (1 — z/|)/£y 
in contact with a rigid body having a surface roughness profile which result in the same 
undeformed gap between the surfaces. 

One of us (Persson) has recently developed a theory of contact mechanics [21], valid 
for randomly rough (e.g., self afnne fractal) surfaces. In the context of rubber friction, 
which motivated this theory, mainly elastic deformation occurs. However, the theory 
can also be applied when both elastic and plastic deformations occur in the contact 
areas. This case is of course relevant to almost all materials other than rubber. 

The basic idea behind the new contact theory is that it is very important not to 
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Figure 20. Three models of "rough" surfaces. In case (a) all the "asperities" are 
equally high and have identical radius of curvature. In this case, according to the 
Hertz contact theory, the area of real contact AA between a solid with a flat surface 
and the surface depends non-linearly on the squeezing force (or load) Fn according to 
'. If the asperities have instead random height distribution as in (b) then, 
for small Fn, AA is nearly proportional to the squeezing force. If the surface roughness 
is random with "asperities" of different heights and curvature radii as in (c), the area 
of real contact for small Fn is exactly proportional to the squeezing force. 

exclude a priori any roughness length scale from the analysis. Thus, if A(X) is the 
(apparent) area of contact on the length scale || A (see Fig. 22), then we study the 
function P(() = A(X)/A(L) which is the relative fraction of the surface area where 
contact occurs on the length scale X = L/C, (where C, > 1), with P(l) = 1. Here 
A(L) = A denotes the macroscopic contact area [L is the diameter of the macroscopic 
contact area so that A ps L 2 ] . 

Consider the system at the length scale X = L/C,, where L is the diameter of the 
nominal contact area. We define q L = 2n/L and write q = q L (. Let P{a,C/) denote the 
stress distribution in the contact areas under the magnification (. The function P(a, () 

| We define the apparent contact area A(X) on the length scale A to be the area of real contact if the 
surface would be smooth on all length scales shorter than A. That is, considering the Fourier transform 
of the surface profile, all the components whose wave vector is larger that 2tt/\ have to be set to 0, 
and the contact area with this new surface is by definition A(X). 




Figure 21. A rubber block (dotted area) in adhesive contact with a hard rough 
substrate (dashed area). The substrate has roughness on many different length scales 
and the rubber makes partial contact with the substrate on all length scales. When a 
contact area is studied at low magnification (£ = 1) it appears as if complete contact 
occurs in the macro asperity contact regions, but when the magnification is increased 
it is observed that in reality only partial contact occurs. 



satisfies the differential equation (see Ref. [21]): 
dP d 2 P 

where f(() = G'(()<j 2 , a being the average pressure in the nominal contact area and 

G(CH^) 2 f'd^Cfa), (6) 

with E* = E/(l - v 2 ). 

Eq. (5) is a diffusion-type equation, where time is replaced by magnification (, and 
the spatial coordinate with the stress a (and where the "diffusion constant" depends on 
(). Hence, when we study P(a,() on shorter and shorter length scales (corresponding 
to increasing (), the P(a,() function will become broader and broader in u-space. We 
can take into account that detachment will occur when the local stress reaches a — 
(we assume no adhesion) via the boundary condition [22]: 

^(0,0 = 0. (7) 

In order to solve the equation (5) we also need an "initial" condition. This is determined 
by the pressure distribution at the lowest magnification ( = 1. If we assume a constant 
pressure a in the nominal contact area, then P(cr, 1) = 5(a — a ). 

We assume that only elastic deformation occurs (i.e., the yield stress a Y — > oo). In 
this case 

POO 

P(C)= / daP(a,C) 
Jo 
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C=L/A, 




P(Q = AW/A(L) 



Figure 22. An elastic ball squeezed against a hard, rough, substrate. Left: the system 
at two different magnifications. Right: The area of contact A(\) on the length scale 
A is defined as the area of real contact when the surface roughness on length scales 
shorter than A has been removed. 



When adhesion is taken into account, tensile stresses can occur at the interface between 
the two solids, and the boundary condition (7) is no longer valid [23], see Sec. 7.1. It is 
straightforward to solve (5) with the boundary conditions P(0, () = and P(oo, () = 
to get 

p(o = 2 r dx !^e- ao (o =erf / r 1 V (8 ) 

V ; n Jo x \2VGJ 

Note that for small load a , G 3> 1 and in this case (8) reduces to P(Q « -Pi(C) where 
P 1 (C) = [vrG(C)]- 1/2 . (9) 

Since G ~ 1 jo\ it follows that the area of real contact is proportional to the load for 
small load. Using (8) and (9) we can write in a general case 



P(C) = erf (^pMC)J (10) 

The physical meaning of the diffusion-like Eq. (5) is as follows: When the system is 
studied at the lowest magnification ( = 1 no surface roughness can be observed and the 
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block makes (apparent) contact with the substrate everywhere in the nominal contact 
area. In this case, if we neglect friction at the interface, the stress at the interface 
will everywhere equal the applied stress <7o, see Fig. 23(a), so that the distribution will 
initially be delta function- like, P(o~,l) = 5(o~ — a ). Increasing the magnification, we 
include surface roughness with wavelength down to A = L/(, and here one may observe 
some non-contact regions as shown in Fig. 23(b). Since the stress must go continuously 
to zero at the edges of the boundary between the contact and non-contact regions, it 
follows that the stress distribution P(a, () will have a tail extending the whole way 
down to the zero stress as indicated in Fig. 23(b) (right). There will also be a tail 
toward larger stresses a > a because the average stress must be equal to uq. This 
distribution broadens as in a diffusion problem. With increasing magnification, the 
stress distribution will broaden further and without limit as indicated in Fig. 23 (right). 

The theory presented above predicts that the area of contact increases linearly with 
the load for small load. In the standard theory of Greenwood and Williamson [18] this 
result holds only approximately and a comparison of the prediction of their theory with 
the present theory is therefore difficult. Bush et al. [20] have developed a more general 
and accurate contact theory. They assumed that the rough surface consists of a mean 
plane with hills and valleys randomly distributed on it. The summits of these hills are 
approximated by paraboloids, whose distribution of heights and principal curvatures 
are obtained from the random process theory. This is to be compared with the GW 
assumption that the caps of the asperities are spherical, each having the same mean 
radius of curvature. As a result of the more random nature of the surface, Bush et al 
found that at small load the area of contact depends linearly on the load accordingly to 



where F N is the normal load, E* = E/(l — u 2 ), and k = (2-ir) 1 / 2 . This result is very 
similar to the prediction of the present theory where, for small load, from (6) and (9), 
A is again given by Eq. (11) but now with k = (8/vr) 1 / 2 . Thus our contact area is a 
factor of 2/n smaller than the one predicted by the theory of Bush et al. Both the 
theory of Greenwood and Williamson and that of Bush et al., assume that the asperity 
contact regions are independent. However, as discussed in Ref. [22], for real surfaces 
(which always have surface roughness on many different length scales) this will never be 
the case even at a very low nominal contact pressure. As we argued [22], this may be 
the origin of the 2/7r-difference between our theory (which assumes roughness on many 
different length scales) and the result of Bush et al. 

The predictions of the contact theories of Bush et al. [20] and Persson [21] were 
compared to numerical calculations (see Ref. [22] [24]). Borri-Brunetto et al. [25] studied 
the contact between self affine fractal surfaces using an essentially exact numerical 
method. They found that the contact area is proportional to the squeezing force for small 
squeezing forces. Furthermore, it was found that the slope a(C) of the line A = a(()F 
decreased with increasing magnification (. This is also predicted by the analytical 
theory [Eq. (11)]. In fact, good agreement was found between theory and computer 




(11) 




Figure 23. The stress distribution P(<r, () in the contact region between a (rigid) 
block and an elastic substrate at increasing magnification (. At the lowest (engineering) 
magnification £ = 1 the substrate surface looks smooth and the block makes (apparent) 
contact with the substrate in the whole nominal contact area. As the magnification 
increases, we observe that the area of (apparent) contact decreases, while the stress 
distribution becomes broader and broader. 

simulations for the change in the slope with magnification and its dependence on the 
fractal dimension Df. 

Hyun et al. performed a finite-element analysis of contact between elastic self-afline 
surfaces. The simulations were done for a rough elastic surface contacting a perfectly 
rigid flat surface. The elastic solid was discretized into blocks and the surface nodes 
form a square grid. The contact algorithm identified all nodes on the top surface that 
attempt to penetrate the flat bottom surface. The total contact area A was obtained by 
multiplying the number of penetrating nodes by the area of each square associated with 
each node. As long as the squeezing force was so small that the contact area remained 
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Figure 24. Dots: the factor function of Hurst's exponent H for self affine 

fractal surfaces. The two horizontal lines are the predictions of the theories of Bush et 
al. (solid line) and Persson (dashed line). From Ref. [24]. 



below 10% of the nominal contact area, A/Aq < 0.1, the area of real contact was found 
to be proportional to the squeezing force in accordance with Eq. (11). Fig. 24 shows 
Hyun et al.'s results for the factor k in (11) as a function of Hurst's exponent H for 
self affine fractal surfaces. The two horizontal lines are predictions of the theories of 
Bush et al. (solid line) and Persson (dashed line). The agreement with the analytical 
predictions is quite good considering the ambiguities in the discretization of the surface. 
The algorithm only consider nodal heights and assumes that contact of a node implies 
contact over the entire corresponding square. This procedure would be accurate if the 
spacing between nodes where much smaller than the typical size of asperity contacts. 
However, the majority of the contact area consists of clusters containing only one or a 
few nodes. Since the number of large clusters grows as H — > 1, this may explain why 
the numerical results approach Persson's prediction in this limit. 

Hyun et al. also studied the distribution of connected contact regions and the 
contact morphology. In addition, the interfacial stress distribution was considered, and 
it was found that the stress distribution remained non-zero as the stress a — > 0. This 
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violates the boundary condition (7) that P(cr, () = for a = 0. However, it has been 
shown analytically [22] that for "smooth" surface roughness this latter condition must be 
satisfied, and we believe that the violation of this boundary condition in the numerical 
simulations may reflect the way the solid was discretized and the way the contact area 
was defined in Hyun et al.'s numerical procedure. 

Elastic contact theory and numerical simulations show that in the region where the 
contact area is proportional to the squeezing force, the stress distribution at the interface 
is independent of the squeezing force. In addition, for an infinite system the distribution 
of sizes of the contact regions does not depend on the squeezing force (for small squeezing 
forces). Thus, when the squeezing force increases, new contact regions are formed in 
such a way that the distribution of contact regions and the pressure distribution remains 
unchanged. This is the physical origin of Coulomb's friction law which states that the 
friction force is proportional to the normal (or squeezing) force [14], and which usually 
holds accurately as long as the block-substrate adhesional interaction can be neglected 
[26]. 

4-2. Surface stiffness of fractal surfaces 

The contact behaviour of realistic surfaces with random multiscale roughness remains 
largely unknown. Recently experimental results of the surface stiffness for self-affine 
fractal surfaces were presented [27]. A hard micrometric flat probe surface was squeezed 
against different substrates, and the vertical displacement u was measured as a function 
of the squeezing force F. It was found that the surface stiffness S = alF/alu depends 
remarkably on the fractal dimension of the substrate surface, decreasing by more than 
an order of magnitude as the substrate fractal dimension is increased by only 10%. 
As expected, the surface stiffness was nearly proportional to the inverse of the rms- 
roughness amplitude a. Fig. 25 shows the surface stiffness as a function of the squeezing 
force F for fractal substrate surfaces produced from two different materials (+ and x 
symbols, respectively). Note the strong decrease in the surface stiffness S with increasing 
fractal dimension _D f of the substrate. 

4-3. Viscoelastic contact mechanics 

The contact between a viscoelastic solid and hard, randomly rough, substrates is a topic 
of great practical importance, e.g., for pressure sensitive adhesives, rubber friction and 
rubber seals. When a viscoelastic solid is squeezed with a constant force against a rough 
substrate, the area of real contact will increase monotonically with the contact time, 
see Fig. 26. Since rubber-like materials have a wide distribution of relaxation times, the 
area of real contact will usually increase over a very long time period (which, e.g., could 
be a year or more). Since the pull-off force depend on the area of real contact, contact 
theories for viscoelastic solids are important for estimating how the pull-off force (or 
tack) depends on the applied squeezing pressure and the squeezing time, see Sec. 4.4. 
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Figure 25. Surface stiffness 5 = dF/du as a function of the squeezing force F for a 
flat hard micrometer sized surface in contact with self affine fractal substrates made 
from two different materials (+ and x symbols, respectively). The fractal dimension 
Df and the rms-roughness amplitude a of the surfaces are: (a) Df = 2.10, a = 20 nm; 
(b) 2.30, 80 nm; (c) 2.15, 36 nm and (d) 2.26, 160 nm. From Ref. [27] 



Rubber-like materials have elastic modulus E(u) which depend strongly on 
frequency uj and temperature T. Thus, at very low frequencies or high temperatures 
they behave as very soft rubbery materials, with typical elastic modulus in the range 
0.01 — 1 MPa. At high frequencies or low temperatures they instead behave as hard 
glassy materials with the elastic modulus of order 1 GPa or more. Thus, as a function 
of increasing frequency (or decreasing temperature) the elastic modulus may increase 
by a factor of 1000 or more. The transition from the rubbery region to the glassy region 
is very wide, usually extending over more than three frequency-decades. In a contact 
experiment, the inverse of the contact time is a characteristic frequency; thus for long 
contact times rubber behaves as a soft solid and the contact area is large, while for short 
contact times it is relatively hard and the contact area is small. 

When a viscoelastic solid is squeezed with a constant force against a hard rough 
substrate, the contact area will increase continuously with time [28, 29]. The relative 
contact area P((,t) = A((,t)/A is approximately given by Eq. (10): 



(12) 



7T 



P(C,t) = erf P- p i(C,i) 





Figure 26. With increasing time, the contact area between the viscoelastic substrate 
and a flat hard solid surface increases through growth of existing asperity contact areas 
(A and C), as well as by the formation of new asperity contact areas (B). 
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where [28] 

p ^-^mL^W) e ~" (13) 

denotes the relative contact area to linear order in <7q. Here &o(uj) is the Fourier 
transform of the squeezing pressure a (t), and 

Q(C) = ( [J dg q 3 C(q)\ . 
Now, let us assume that u (t) = o"o for < t < t\ and zero otherwise. In this case 

^L(l-e^) (14) 
Substituting this into (13) gives 
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Figure 27. The relative contact area [obtained from Eqs. (12) and (15)] at the 
magnification £ = 10 5 , as a function of the logarithm of the contact time t, for the 
nominal pressures cro = 0.1 MPa, for a tread tire rubber compound, and for a compound 
used in the tire-rim area, at the temperature T = 60° C. The substrate is self affinc 
fractal with the fractal dimension Df = 2.2, and with the rms roughness 1 /zm, and with 
the roll-off wave vector (see Fig. 6) go = 6 x 10 4 m _1 (this is a typical qo for a polished 
steel surface). The magnification ( = 1 refer to the length scale Ao = 2n/qo ss 0.1 mm, 
so that the magnification ( = 10 5 correspond to the length scale (resolution) A w 1 nm. 



This result does, of course, not depend on t\ as long as t < t\ (causality). Thus we are 
free to choose for t\ any value larger than the time t under consideration. The equations 
presented above are only valid as long as the contact area increases with time, which is 
the case, e.g., if a constant squeezing force is applied at time t — 0. As an illustration, 
in Fig. 27 we show the relative contact area calculated from Eq. (15) using the measured 
viscoelastic modulus for two types of rubber. The rubber block is squeezed with the 
nominal pressure o"o = 0.1 MPa against a "steel" surface. We show results for two 
typical rubbers used in the construction of tires, namely a tread rubber compound, and 
a rubber compound from the tire rim region involved in the tire air sealing (see Sec. 5). 
The substrate was a self affine fractal with the fractal dimension Df = 2.2, and with the 
rms-roughness amplitude 1 /xm. The tread compound exhibits the largest contact area 
because it is elastically softer than the tire rim compound. Note also that the change in 
the contact area is only a factor ~ 2 — 4 when the contact time changes by ~ 9 decades. 




(15) 
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44. Tack 

Pressure sensitive adhesives (PSA) are used in many important applications, e.g., for 
Scotch tapes, post-it pads and self-adhesive labels and envelopes [30, 31]. The adhesive 
consists of a very thin layer (usually of order ~ 20 — 100 /im) of a very soft, weakly 
cross-linked rubber compound. The low frequency elastic modulus is typically only 
0.01 MPa which is ~ 100 times lower than the rubber used for tires. As a result of 
the low elastic modulus, nearly complete contact will occur in the apparent contact 
area even for relatively low squeezing pressures and large surface roughness. This is in 
contrast to tire rubber which under similar condition would give a contact area of only 
a few % of the nominal contact 3X6 cl ? 3S discussed in Sec. 6.1. 

We carried out calculations [using Eq. (15)] of the time-dependent contact area for 
a standard rubber tack compound [acrylic PSA namely poly(2ethylhexyl acrylate) with 
2% acrylic acid (PEHA-AA)]. 

Fig. 28 shows the results at T = 20° C when the PSA is squeezed in contact with 
"steel" surfaces with rms roughness amplitudes ho = 0.25, 1 and 4 /im. Fig. 28(a) shows 
the logarithm of the relative contact 3r<33. 3S 3 function of the logarithm of the contact 
time, for a contact pressure o"o = 0.01 MPa. Fig. 28(b) shows the logarithm of the 
relative contact area after Is of contact, as a function of the logarithm of the contact 
pressure. 

It is interesting to compare the results in Fig. 28 with the experimental data 
reported in Ref. [32] and shown in Fig. 29. The data show the dependence of the 
maximum pull-off force F on the contact time and contact pressure for smooth and 
rough PSA on the same smooth steel surface. The tack film is a standard polymer 
compound similar to PEHA. In Ref. [32] no numerical values of the rms roughness was 
presented for any of the studied surfaces. If one assumes, as a first approximation, 
that the pull-off force is proportional to the area of real contact, and moreover that the 
smooth and rough PSA films in the measurements correspond, say, to the surfaces in 
Fig. 28 with the rms roughness 0.25, 1 and 4 /im, then the agreement between theory 
and the experiment is remarkably good. 

5. Seals 

Surface roughness is an important factor which influences the rate of leakage through 
seals. The exact mechanism of roughness induced leakage is not well understood [33]. 
In this section we present a new way of looking at this problem [28]. 

Viscoelastic materials such as rubber are often used for sealing. Here we consider 
the tire-rim sealing. We are interested in the air (or gas) flow from inside the tire to the 
outside via the roughness-induced non-contact area (pore channels) in the rubber-steel 
rim area. The rim is made from steel. We assume that the steel surface is a self-afnne 
fractal with the fractal dimension Df = 2.2. The surface root-mean-square roughness is 
assumed to be 1 /im, a rather typical value for polished steel surfaces. 
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Figure 28. (a) Calculated logarithm of the relative contact area at the magnification 
C = 10 5 , as a function of the logarithm of the contact time. Contact pressure 
(To = 0.01 MPa. (b) Calculated logarithm of the relative contact area after Is of 
contact at the magnification ( = 10 5 , as a function of the logarithm of the contact 
pressure. All calculations for PEHA-AA at temperatures T = 20° C on a self affinc 
fractal substrate with fractal dimension Df = 2.2, roll-off wave vector q n = 6 x 10 4 m _1 
and rms roughness amplitudes ho — 0.25, 1 and 4 fim. 



Fig. 30 shows the calculated [using Eq. (15)] relative contact area after one second 
of contact for tread tire rubber (a), and for rim tire rubber (b), and the root-mean- 
square roughness (c), as a function of the magnification (the time dependence of the 
contact area was presented in Fig. 27). Here we defined 

/ rqi \ V 2 

/W(C)= 27r / d <7<?C(<?) 

The results are for a nominal pressures o-q = 0.1, 0.2 and 0.4 MPa, at the temperature 
T = 60°C. 

We consider now a rubber-steel interface at increasing magnification. At the lowest 
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Figure 29. (a) Logarithm of the experimental pull-off force F (in arbitrary units) 
as a function of the logarithm of the contact time. Contact pressure Co = 0-01 MPa. 
From Ref. [32] (b) Logarithm of the experimental pull-off force [same units as in (a)] 
after Is of contact, as a function of the logarithm of the contact pressure, for three 
different surface roughnesses. Ref. [32] 



magnification ( = 1 complete contact occurs at the interface, see Fig. 31 (left). When 
we increase the magnification we observe non-contact areas or islands. Magnification 
is now increased until the non-contact areas percolate, i.e., until a channel of non- 
contact surface area, extending from the high-pressure internal region of the tire to 
the outside (atmospheric pressure region), is first observed. As the magnification is 
increased further, new non-contact region will be observed, but the separation between 
the surfaces in these new non-contact areas will be smaller than along the percolation 
channel. Since the gas flow N (number of molecules per unit time) through a rectangular 
pore of height h depends as iV ~ h 3 we will assume that most gas will leak through the 
percolation channel. 

Assuming that the contact regions at any magnification are approximately 
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Figure 30. Calculated relative contact area after one second of contact for tread tire 
rubber (a), for rim tire rubber (b), and the root-mean-square roughness as a function 
of magnification (c), where ( = 1 corresponds to a wavelength A = 2ir/q ss 100 ^m. 
The highest magnification ( = 10 5 corresponds to A w 1 nm. Results are shown for a 
nominal pressures ctq = 0.1, 0.2 and 0.4 MPa, at the temperature T = 60°C. Substrate 
surface with rms roughness ho = 1 /im and roll-off wave vector qo = 6 x 10 m _1 . 
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Figure 31. When the interface between the solids is studied at low magnification, 
there appear to be complete contact (black area) between the solids. When the 
magnification is increased it is observed that only partial contact occur. At high 
enough magnification the non-contact (white) surface area will percolate and one (or 
several) air flow channels will be visible at the interface. 



randomly distributed in the apparent contact region, we expect from percolation theory 
that the non-contact region will percolate when A/A ~ 1 — p c , where p c is the site 
percolation number [34]. For a hexagonal lattice, which is the most plausible lattice 
structure in the present case, one has p c ~ 0.7, while for a square lattice (as in Fig. 31) 
percolation occur at p c ~ 0.6. Thus, the exact value of the percolation threshold does 
not depend sensitively on the symmetry of the unit cell. 

We assume that the main gas leakage comes from gas flow through the percolation 
channel. The narrowest passage in this channel can be considered as a rectangular pore 
of height h, and of width and length A, where A is determined by the magnification ( c 
at the point where A/A = 1 — p c pj 0.3. The height h of the pore is (approximately) 
determined by the rms roughness at the magnification ( c . In the present case, if the tire 
gas pressure is in the range 0.2 — 0.3 MPa, from Fig. 30(b) we get ( c ~ 10 and from 
Fig. 30(c), h « 0.1 fim. 

We divide the tire-rim contact area into m square areas B x B, where B is the 
width of the tire-rim contact area (we expect B to be of order a few cm). The number 
of squares is m = 2nR/B, where R is the radius of the tire at the rim. We expect 
m m 100. 

Let us study the gas flow through a rectangular pore of height h and width and 
length (in the flow direction) A. We assume stationary and laminar flow, and that 
/»< A. In this case, if Ni(t) denote the number of gas phase molecules in the tire (which 
is proportional to the pressure Pi in the tire), the basic equations of hydrodynamics give 
for the typical case P 1 3> Pq- 

mP?h 3 

iVi ~ 

24fik B T 

Here we have implicitly assumed that the full pressure drop P\ — Pq occurs over the 
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pore. Thus, the time it takes for the pressure in the tire to drop with 4% is 
At = — — 

where V\ = Vq/itl (where Vq is the air volume in the tire) is the volume of air in an 
angular section of the tire of width B. With B tv 3 cm we get V\ ~ 3 x m 3 and 
using the viscosity of air /j, 17 x 1CT 6 Ns/m 2 gives At m 1 year. This is an upper 
limit of the leakage time, since when the interfacial contact area is studied at higher 
magnification new pore channels through which the air can leak will be detected. 

It is interesting to note that the adhesional interaction between the rubber surface 
and the steel rim is likely to have negligible influence on the leakage rate. Adhesion will 
affect the (apparent) contact area only at very high magnification [28], but most of the 
gas leakage occurs via the much larger air flow channels which can be observed at low 
magnification. 

When a rubber block is squeezed with a constant force against a rough substrate, 
the area of real contact will increase continuously, and the (average) space between 
the surfaces will decrease continuously with increasing time. This is due to stress 
relaxation in the rubber and was illustrated in Fig. 27. In some applications, such 
as tire-air sealing, this relaxation effect is clearly beneficial as it will reduce the size of 
the interfacial air-flow channels. However, in other applications, such as O-ring seals, 
stress relaxation may result in catastrophic events. One prominent example was the 
Challenger catastrophe [35]. In this case rubber O-rings where used to seal the hot gas 
inside the solid fuel rockets and prevent it from leaking through the horizontal joints used 
to hold the different vertical rocket wall-segments together. O-rings normally operate 
with about ~ 15% compression to ensure a tight seal, see Fig. 32 (left). However, the 
high pressure inside the rocket resulted in expansion of the steel tube, and to a (small) 
separation of the wall-surfaces at the joints. At high enough rubber temperature this 
is not a problem since the compressed rubber O-ring would then quickly expand and 
seal tightly even if the space between the walls would increase, see Fig. 32 (middle). 
However, the temperature at the time of the Challenger launch was exceptionally low 
(around zero Celsius degrees), which resulted in very long rubber relaxation times, so 
that the compressed form of the O-seals was "frozen in", see Fig. 32 (right). This effect 
was demonstrated by Feynman in a famous experiment where part of an O-ring was 
kept in a deformed state in a glass with ice water. Thus, as the space between the 
surfaces increased, one or several O-rings failed to seal-off the hot gas, resulting in leak 
of hot fuel gas through the rocket joints, and finally to the catastrophic failure of the 
rocket. 

Stress relaxation is an important aspects of O-ring seals even when no separation of 
the solids walls occurs at the seal [36]. When a seal is under constant compression (fixed 
solid walls), the initial stress decays with increasing time, roughly proportional to the 
logarithm of the time of contact. Thus the peak compressive stress may eventually drop 
below the system (gas or liquid) pressure, and the seal leaks. Clearly, stress relaxation 
effects must be taken into account when determining the type of O-ring to be used in a 
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Figure 32. Rubber O-ring seal. Left: The rubber ring is compressed between two 
solid walls. Middle: If the rubber temperature is high enough, when the separation 
between the walls is suddenly increased the rubber will quickly expand and still provide 
tight sealing. Right: At low temperature the rubber relaxation times are very long 
and the rubber will not be able to expand quickly enough. Here the sealing fails. 

particular application. 
6. Rubber friction 

First principle calculations of frictional forces for realistic systems are generally 
impossible. The reason is that friction usually is an interfacial property often 
determined by the last few uncontrolled monolayers of atoms or molecules at the 
interface. An extreme illustration of this is offered by diamond. The friction between two 
clean diamond surfaces in ultra high vacuum is huge because of the strong interaction 
between the surface dangling bonds. However, when the dangling bonds are saturated 
by a hydrogen monolayer (as they generally are in real life conditions), friction becomes 
extremely low [37]. Since most surfaces of practical use are covered by several monolayers 
of contaminant molecules of unknown composition, a quantitative prediction of sliding 
friction coefficients is out of the question. An exception to this may be rubber friction 
on rough surfaces, which is the subject we address here. 

Rubber friction on smooth substrates, e.g., on smooth glass surfaces, has two 
contributions, namely an adhesive (surface) and a hysteretic (bulk) contribution [38, 39]. 
The adhesive contribution results from the attractive binding forces between the rubber 
surface and the substrate. Surface forces are often dominated by weak attractive van 
der Waals interactions. For very smooth substrates, because of the low elastic moduli 
of rubber-like materials, even when the applied squeezing force is very gentle this weak 
attraction may result in a nearly complete contact at the interface [23, 40], leading to 
the large sliding friction force usually observed [41]. For rough surfaces, on the other 
hand, the adhesive contribution to rubber friction will be much smaller because of the 
small contact area. The actual contact area between a tire and the road surface, for 
example, is typically only ~ 1% of the nominal footprint contact area [21, 42, 43]. Under 
these conditions the bulk (hysteretic) friction mechanism is believed to prevail [21, 43]. 
For example, the exquisite sensitivity of tire-road friction to temperature just reflects 
the strong temperature dependence of the viscoelastic bulk properties of rubber. 

Here we discuss how rubber friction depends on the surface roughness power spectra, 
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and we consider the influence of wear-induced polishing and of water on the road track 
on rubber friction. 

6.1. Basic theory of rubber friction 

The main contribution to rubber friction when a rubber block is sliding on a rough 
substrate, such as a tire on a road surface, is due to the viscoelastic energy dissipation 
in the surface region of the rubber as a result of the pulsating forces acting on the rubber 
surface from the substrate asperities, see Fig. 33. Recently one of us has developed a 
theory which accurately describes this energy dissipation process, and which predict 
the velocity dependence (and, in a more general case, the time-history dependence) 
of the rubber friction coefficient [21, 42]. The results depend only on the (complex) 
viscoelastic modulus E(cu) of the rubber, and on the substrate surface roughness power 
spectra C (q) . Neglecting the flash temperature effect (the term flash temperature refers 
to a local and sharp temperature rise occurring in the tire-road asperity contact regions 
during slip), the kinetic friction coefficient at velocity v is determined by [21] 



where a is the mean perpendicular pressure (load divided by the nominal contact area), 
and v the Poisson ratio, which equals 0.5 for rubber-like materials. 

The theory takes into account the substrate roughness in the range q$ < q < q*, 
where q^ is the smallest relevant wave vector of order 2n/L, where (in the case of a 
tire) L is the lateral size of a tread block, and where q\ may have different origins (see 
below). Since ql for a tire tread block is smaller than the roll-off wave vector q of the 
power spectra of most road surfaces (see Fig. 15), rubber friction is very insensitive to 
the exact value of q^. 

The large wave vector cut off q\ may be related to road contamination, or may 
be an intrinsic property of the tire rubber. For example, if the road surface is covered 
by small contamination particles (diameter D) then q\ tv 2tt/D. In this case, the 
physical picture is that when the tire rubber surface is covered by hard particles of 
linear size D, the rubber will not be able to penetrate into surface roughness "cavities" 
with diameter (or wavelength) smaller than D, and such short-range roughness will 
therefore not contribute to the rubber friction. For perfectly clean road surfaces we 
believe instead that the cut-off g* is related to the tire rubber properties. Thus, the 
high local (flash) temperatures during braking, and the high local stresses which occur 
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Figure 33. Hysteretic friction of a rubber block on a rough road surface. The 
road asperities exert pulsating forces on the sliding rubber block, leading to energy 
dissipation in the rubber via the rubber internal friction. Most of the energy dissipation 
occurs in the volume elements bounded by the dashed lines. The rubber viscoelastic 
deformations in the large volume elements are induced by the large road asperities, 
while the smaller dissipative regions result from the smaller asperities distributed on 
top of the large asperities. In general, in calculating the rubber friction, the viscoelastic 
energy dissipation induced by all the asperity sizes must be included, and the local 
temperature increase (flash temperature) in the rubber resulting from the energy 
dissipation should also be taken into account in the analysis. 

in the tire rubber-road asperity contact regions, may result in a thin (typically of order 
a few micrometer) surface layer of rubber with modified properties (a "dead" layer), 
which would contribute very little to the observed rubber friction. Since the stresses 
and temperatures which develop in the asperity contact regions depend somewhat on 
the type of road [via the surface roughness power spectra C(q)} ) the thickness of this 
"dead" layer may vary from one road surface to another, and some run-in time period 
will be necessary for a new "dead" layer to form when a car switches from one road 
surface to another. Such "run-in" effects are well known experimentally. 

In the theory that one of us has developed, the thickness of the dead layer is 
determined by studying (via computer simulations) the temperatures and stresses which 
develop on the surfaces of the tire tread blocks during ABS breaking. 

Since this is not the proper place for a full presentation of all details of the 
theoretical calculations of tire-road friction, we shall simply present some numerical 
results to illustrate how the rubber friction depends on the surface roughness power 
spectra. Fig. 34(a) shows the kinetic friction coefficient as a function of the sliding 
velocity calculated for a rubber block sliding on a self afnne fractal surface with the 
Hurst exponents H = 0.8, roll-off wave vector q = 1000 m -1 , and large wave vector 




Figure 34. Calculated friction coefficients for a rubber block sliding on a self affinc 
fractal surfaces with H = 0.8, roll-off wavevector q = 1000 m -1 , largest wave vector 
Qi = 1o x 320.0 (corresponding to the magnification £ =320) rms roughness a = 1 mm, 
and the background temperature assumed was To = 60°C. (a) The kinetic friction 
coefficient as a function of the logarithm of the sliding velocity. Solid and the 
dashed lines are with and without flash temperature effects, respectively, (b) Flash 
temperature calculated at two different magnifications. The highest magnification 
C = 300 correspond to the temperature very close to the surface while the low 
magnification curve correspond to the temperature deeper into the rubber, (c) Relative 
contact area P(Q = A(()/A calculated for two different magnifications. 
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cut-off q{ = 320.0 x q , typical for road surfaces. We show results both with and without 
the flash temperature effect. 

For low sliding velocities, the kinetic sliding friction depends very weakly on the 
sliding velocity, and only at extremely low velocities is friction strongly reduced. Thus 
the present calculation (as well as experiments) shows that even sliding velocities as small 
as 10~ 10 m/s, may give rise to a relative large kinetic friction. The physical reason for 
this is the very wide distribution of relaxation times in most rubber materials, extending 
to very long times. 

For typical tire rubber compounds, calculations (neglecting temperature effects) as 
well as measurements at low sliding velocities and different temperatures (subsequently 
shifted to a common temperature utilizing the frequency-temperature Williams-Landel- 
Ferry (WLF) transform [44]), indicate that the maximum friction coefficient ^(v) will 
typically occur at very high velocities, of order 1000 m/s [see dashed line in Fig. 34(a)]. 
However, direct experiments without the WLF transform show a maximum in fi^ at 
much lower velocities. This is a result of the flash temperature, see solid line in Fig. 34(a). 
Thus when the flash temperature effect is taken into account the maximum of fJ>k(v) shifts 
at much lower sliding velocities in the typical range 0.1 — 1 cm/s, in agreement with 
experiment. 

Fig. 34(b) and (c) show the flash temperature and the relative contact area at 
two different magnifications. For velocities v < 10~ 4 m/s the contact area decreases 
with increasing sliding velocity. This is a result of the increasing frequencies of 
the deformations of the tire surface with increasing sliding velocity. However, for 
v > 10~ 3 m/s the contact area increases. This is a result of the increased temperature of 
the rubber in the contact areas shifting the viscoelastic dissipation maximum to higher 
frequencies and hence making the rubber elastically softer at any given frequency. It 
is also of interest to consider the relative contact area. For sliding velocities of order 
1 cm/s the tire-road contact area at a magnification of order 100 — 1000 is just a few 
% of the nominal footprint contact area. This is illustrated in Fig. 35 which shows the 
tire-road contact region at different magnification. 

The flash temperature effect described above is of extreme practical importance. If 
the excessive buildup of flash temperatures could be avoided or reduced, the effective 
tire-road friction coefficient would clearly increase. One way to reduce the flash 
temperature is to use wide tires instead of narrow tires^f. This can be understood 
as follows. In a first approximation the nominal tire-road contact area is the same for 

The discovery that wide tires gives higher effective friction than narrow tires was made by Jim Hall 
and coworkers in the middle 60's for racing tires. The story is beautifully described in [45], where Hall 
stated: "In a matter of less than a year we went from those narrow little Dunlop Green Sports to a tire 
that was almost double the width. And our single-speed torque converter just couldn't handle all that 
grip. We immediately went to a two-speed transmission so we could get more torque multiplication 
off the slow-speed corners." "So that's what happened with me on tires. In the middle '60s, probably 
'64 or '65, Firestone took us from less than 6 inches of tread to width up to 12 inches." Paul Haney 
gives in his book [45] an explanation for why wide tires give higher effective friction than narrow tires. 
However, this 'explanation', and many others given on the Internet, are incorrect. 
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Figure 35. A tire in contact with a road surface does not make perfect contact with 
the road in the whole footprint area. At the lowest magnification the asperity contact 
regions have diameters of order ~ 0.3—1 cm. When the magnification is increased these 
macrocontact regions break up into smaller contact regions separated by noncontact 
areas. At a magnification of order <~ 100 — 1000 the tire rubber-road contact area is 
just a few % of the nominal footprint contact area. 

narrow and wide tires. In fact, this holds exactly when a tire can be described as a thin 
elastic (torus shaped) membrane (which is a good approximation for airplane tires but 
less so for car tires) and the tire air pressure and the load are the same for both tires" 1 ". 
Thus, wide tires will have a footprint area which is shorter in the longitudinal rolling 
direction than a narrow tire, see Fig. 36. In order to fully build up the flash temperature, 
a tread block must slide a distance of order the average diameter D of the macro asperity 
contact regions. The macro asperity contact regions are the contact between the rubber 
and the largest road asperities involving the longest wavelength roughness components 
characterized by the roll-off wavelength A = 27r/g in the surface roughness power 

+ For the membrane tire the nominal pressure in the footprint contact area is given by the air pressure 
P in the tire so that the nominal contact area Aq — Fn /P is the same for both tires if the load Fn is 
the same. 
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Figure 36. The nominal tire-road contact area is approximately the same for a wide 
and a narrow tire. Hence the length w of the foot print area in the rolling direction 
will be smaller for the wide tire. 

spectra (see Fig. 15). In a typical case A ~ 1cm, and the average diameter of the 
macro asperity contact regions D 0.5 cm. If a tread block in the footprint contact 
region slips by less than the distance D, the local temperature increase will be below 
that expected during stationary sliding at the same slip velocity. If v denotes the average 
slip velocity of a tread block then the slip distance d ~ vt where the time the tread 
block spend in the footprint area t = w/v R , i>r being the tire rolling velocity, and w the 
length of the footprint area. Thus, d « (v/vr)w and the slip velocity necessary for the 
slip to be of order D is v (D/w)v R . Hence, for a wide tire (small w) the tread block 
slip velocity can reach higher values than for a narrow tire without building up the flash 
temperature effect fully. Since the kinetic friction coefficient in the absence of the flash 
temperature effect increases monotonically with increasing slip velocities (up to very 
high slip velocities, perhaps 1000 m/s or more, see dashed line in Fig. 34) it follows that 
wide tires will give higher effective friction than narrow tires. The discussion above is 
only of semi- quantitative validity since in reality the tread blocks will not slip with a 
constant velocity in the footprint area, but with a nonuniform velocity. 

Another way to reduce the flash temperature could be to use a filler material which 
exhibits a first order phase transition somewhere close to (but above) the background 
tire-rubber temperature, which typically may be around 60°C. If the phase transition 
(e.g., a structural, magnetic or ferroelectric transition) could absorb a sufficiently large 
latent heat, it would strongly reduce the flash temperature. The filler materials used 
today (mainly carbon black and silica) does not exhibit any such phase transition. 

The theoretical results presented above are in accordance with experimental 
observations [46, 47, 48, 49]. Thus, flash temperature effects have been observed during 
dry ABS braking test, where the temperature distribution on the tire surface can be 



On the nature of surface roughness with applications 



43 



measured with high-speed IR-cameras [46, 50]. In one experiment the temperature 
increase in the hot spots on the tread block surfaces on the length scale of the 
measurement was ~ 25 C above the background temperature. The decrease of the 
kinetic friction force with increasing sliding velocity for v > 0.1 — lcm/s (as a result 
of the flash temperature), has been reported by several research groups [46, 47, 48, 49] 
and is in good agreement with the theory. A decreasing kinetic friction coefficient may 
result in stick-slip motion; thus one mechanism of rubber stick-slip is due to the flash 
temperature. 

The theory also predicts that when a rubber block is sliding on a rough substrate 
the friction force is nearly proportional to the normal load, i.e., the friction coefficient is 
independent of the load. This is again in agreement with experiment, see, e.g., Ref. [47]. 
We note that for tires the effective friction coefficient usually decreases with increasing 
load. This does not actually reflect a load dependence of the fundamental rubber- 
substrate friction coefficient, but results rather as an effect of the dependence of the size 
and the pressure distribution in the footprint area on the load, which affects the motion 
of the tread blocks [51]. 

The theory sketched above has been used successfully to predict the ABS breaking 
performance of tires [52, 46, 50, 49]. In one study 6 tires were produced with variation 
only in the tread compound [46, 50]. The tread rubber compounds were varied so that 
big variations in ABS rating were expected. The observed braking distance was found 
to correlate very well with the theoretical predicted effective friction coefficient, with a 
high correlation coefficient R 2 = 0.97. This result shows that it is possible to estimate 
the performance of a tire on physically sound frictional models which focus not only 
on the viscoelastic properties of the rubber compound, but take explicitly into account 
additional important information like the surface topography at all length scales as well 
as ambient and flash temperature conditions. 

6.2. Rubber friction and the influence of polishing 

Let us now discuss the role of road polishing by sliding tires*. It is known that a relative 
new asphalt roadway has a "sharper" surface pavement with higher coefficient of friction 
compared to an older well-traveled asphalt roadway. This is a consequence of the wear 
action exerted by vehicle tires, which results in polishing the exposing aggregate surface, 
and the state of polish is one of the main factors affecting the resistance to skidding. 
(Note: the fact that the polishing of the surface asperities reduces the friction is in 
itself an indication that the adhesional contribution to rubber friction on road surfaces 
is unimportant, since the latter should increase when the surface becomes smoother.) 
Resistance to this polishing action is determined principally by the inherent quality of 
the aggregate itself. Rocks composed of minerals of widely different hardness, and rocks 
that wear by the pulling out of mineral grains from a relative soft matrix, have high 

* See, e.g., http://www.betterroads.com/articles/oct03b.htm for experimental information from road 
engineers about rubber friction and road surfaces, in particular about road polishing. 
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resistance to polishing, see Fig. 37(c)-(e). Conversely, rocks consisting of minerals having 
nearly the same hardness wear uniformly and tend to have lower resistance to polishing, 
see Fig. 37(a), (b). Thus, sandstones have a good resistance to polishing, whereas 
the limestone and flint groups present the lowest resistance^. Among other groups, 
basalt, granite and quartzite yield intermediate results. Microroughness is gradually 
polished away by the action of heavy traffic so, over time, the skidding resistance of 
a road will fall to an "equilibrium" level that depends upon the type of aggregates 
used in the surfacing. For concrete and asphalt pavements that have too low friction, 
diamond grinding will bring its skid resistance up. Grinding increases the pavement 
microstructure by dislodging polished sand particles in the mortar matrix. 

Experience by road engineers shows that rubber friction on dry clean roads depends 
on micro roughness with wavelength mainly in the range 1 fim < A < 1000 //m, in 
accordance with our calculations. Longer wavelength roughness 0.1cm < A < 5 cm is 
important on wet road surfaces especially to prevent hydroplaning for velocities above 
60km/h. High dry friction on asphalt or concrete road surfaces are obtained when 
hard minerals (such as feldspat and quartz) and aggregates with good microstructure, 
such as sandstone and slag (usually, air-cooled blast furnace slag from iron production), 
are embedded in a softer matrix, e.g., tar (asphalt) or mortar (concrete). Minerals 
with rough grains, or mixture of minerals with different texture, will resist polishing 
and maintain high frictional properties. Special surface coatings with very high skid 
resistance have been developed. One such coating consist of epoxy adhesive binding 
a synthetic aggregate to the road surfacesff. The aggregate, manufactured from steel 
slag, has extremely high polish resistance, and is claimed to result in tire-road friction 
about 40% higher than on normal road surfaces. In temperate climates such as Sweden, 
the surface microroughness varies cyclic, with skidding resistance at its lowest during 
the summer and autumn, recovering to some extent during the winter. This is due in 
part to weather induced erosion of the road surface (e.g., freezing water in the road 
cavities and frost heave may break-off road surface fragments), and partly due to the 
strong wear-induced surface roughening which result from the "heavy" equipments used 
to clean the road surfaces from snow and ice. 

Summarizing, it is clear that the best road pavements, with high resistance against 
polishing, consist of aggregates of very small (micrometer) and hard particles, e.g., 

jt The influence of polishing of road surfaces is studied in a standardized lab test where fragments 
(of diameter ~ 1 cm) of the stone (or other material) to be tested are glued onto a flat surface and 
exposed to accelerated (as compared to polishing on the road) polishing using standardized conditions. 
The rubber friction coefficients obtained after a fixed time period of polishing is then measured. The 
result will depend on the rubber used in the friction test (and the laboratory conditions, e.g., the 
temperature), and only relative friction values (for different polished surfaces) are analyzed. In one 
such set of measurements (see http://www.wainwright.co.uk/technical.htm) it was found that gritstone 
and sandstone have the best resistance against polishing (giving friction coefficients after polishing 
/x = 0.68 and 0.67, respectively) followed by basalt and granite (fi = 0.56 and 0.55 respectively), while 
limestone gave the worst result (fi = 0.37). 

ff See http://www.new-technologies.org/ECT/Civil/italgrip.htm for high friction coatings for road 
surfaces. 





Figure 37. (a) and (b): Rolling and braking on cobble stone road surfaces causes the 
surface asperities to get polished, (c)-(e) On asphalt or concrete road surfaces smaller 
hard stone particles never get highly polished as they are continuously exposed and/or 
removed by wear of the softer surrounding matrix. 



fragments of sandstones, where the wear occurs by removal of the individual micro- 
particles or clusters of micro-particles, rather than by polishing the particles, see Fig. 38. 
However, the binding of the particles in the aggregates must be so strong that the wear 
of the road pavement is slow, i.e., the wear rate by removal of particles should be only 
slightly larger than the wear rate by polishing of the particles. The optimal case is when 
the wear process results in particle aggregates of various sizes being removed from the 
road surface so that the road surface remains rough at all the length scales above the 
size of the smallest particles. 

6. 3. Rubber friction on wet road surfaces 

For rubber friction on wet rough substrates at low sliding velocities it is known that the 
friction typically drops by as much as 20 — 30% relative to the corresponding dry case 
[53, 54]. Owing to the small contact area, this cannot be the result of a water-induced 
change of adhesion. On the other hand, as will be discussed below, the friction decrease 
cannot be blamed on a purely hydrodynamical effect either. That leaves finally the 
possibility that water might change precisely the bulk, hysteretic friction. We proposed 
recently [55] that this is indeed the case. Water pools that form in the wet rough 
substrate are sealed off by the rubber, as sketched in Fig. 39, and that will effectively 
smoothen the substrate surface. Smoothening reduces the viscoelastic deformation from 
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Figure 38. Asphalt or concrete roads with aggregates of small (micrometer) particles, 
e.g., fragments of sandstone or slag, have a high resistance against polishing if the tire- 
road interaction mainly results in a slow removal of the microparticles in the aggregates, 
preserving an aggregate particle surface that is (on the microscale) constantly rough, 
(a) and (b) show an aggregate before and after wear (schematic). 

the surface asperities, and thus reduces rubber friction. 




Figure 39. A rubber block sliding on a rough hard substrate, (a) On a dry substrate 
the rubber penetrates a large valley and explores the short wavelength roughness 
within. The pulsating rubber deformations induced by the short-wavelength roughness 
contributes to the friction force, (b) On a wet substrate the valley turns into a water 
pool. Sealing of the pool now prevents the rubber from entering the valley. By removing 
the valley contribution to the frictional force, this sealing effect of rubber reduces the 
overall sliding friction. 
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As discussed before, rubber friction from the viscoelastic deformation by the 
substrate asperities is determined by the complex frequency-dependent bulk viscoelastic 
modulus E(oj) of rubber and by the substrate surface roughness power spectrum C(q). 




Figure 40. (a) Optically observed height profile of a dry asphalt road (1.5 cm x 1.5 cm 
area), darker areas corresponding to deeper regions, (b) Calculated wet profile for the 
same area, with water pools (light blue). 

The upper curve in Fig. 41 shows the power spectrum extracted via Eq. (1) from 
the measured height profile h(x). The log-log scale shows that for q > 1600m -1 , 
C(q) drops as a power law, as expected for a self-affine fractal surface. The fractal 
dimension of this surface is determined by the slope of the curve in Fig. 41 and is about 
D{ = 2.2. The root-mean-square roughness can be obtained directly from the height 
profile, h Tms ~ 0.3 mm. 

Consider now a tire rolling and sliding on a wet road surface. At low velocities 
(say v < 60km/h) there will be negligible hydrodynamic water buildup between the 
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tire and the road surface [55]. In essence, if v < (cr/p) 1 / 2 , where a is the perpendicular 
stress in the tire-road contact area and p the water mass density, there is sufficient time 
for the water to be squeezed out of the contact regions between the tire and the road 
surface, except for water trapped in road cavities. The water pools will be sealed off by 
the road-rubber contact at the upper boundaries of the cavities (see Fig. 39). Thus, we 
can focus on the smoothening effect on the road profile caused by the sealing effect. 

Starting from a dry substrate profile h(x) we can numerically build a new wet 
surface height profile h'(x) as shown in Fig. 40(b). The algorithm assumes every valley 
to be filled with water up to the maximum level where the water will remain confined, 
i.e., up to the lowest point of the edge surrounding the pool. Any extra water added to 
the profile of Fig. 40(b) will flow straight out of the square area. Once the size of the 
area considered is at least as large as Ao this construction becomes unique, with no free 
parameter. 

From the water-smoothened height profile h'(x) we obtain a modified power 
spectrum C'(q) shown by the lower curve in Fig. 41. While the fractal power-law 
decay and the roll-off wave vector are essentially the same as for the dry surface, the 
reduction in the power spectrum reflect the effective water-induced smoothening of the 
rough substrate. 
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Figure 41. Surface roughness power spectra C(q) (above) extracted from the 
measured height profile for a dry asphalt road surface, and (below) calculated assuming 
sealing of all pools in the same surface when wet, as in Fig. 40. Note the logarithmic 
scales. 

The sealed-off water in the pools (see Fig. 39) removes the contact with the interior 
of the valley, which smoothen the effective substrate roughness profile. Our basic 
assumption is therefore that when rubber slides on the wet rough surface, the friction 
force will be determined by the modified power spectrum C\q) < C(q). 
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Let us now examine, based on this model, numerical results about tire friction 
on dry and wet substrates, calculated using the hysteretic friction theory presented in 
Ref. [21]. The hysteretic friction coefficient at velocity v is determined by knowledge of 
the rubber viscoelastic modulus E{u) and of the surface roughness spectrum C(q). 

We present results for the friction of a standard tread compound, sliding on the 
asphalt road just characterized. We used the measured rubber complex viscoelastic 
modulus (not shown) along with the power spectra presented in Fig. 41 for the dry and 
for the wet road surfaces. 



1.2 
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Figure 42. Kinetic friction coefficient as a function of the logarithm of the sliding 
velocity, calculated for a standard tread compound and an asphalt substrate with the 
roughness spectra of Fig. 41. 

Fig. 42 shows the rubber-asphalt kinetic friction coefficient calculated for the dry 
surface (here including also the flash temperature) at T = 60° C as a typical tire 
temperature while rolling on a dry road, and for the wet surface at four different 
temperatures, namely T = 30, 40, 50, and 60°C (on a wet road the tire temperature is 
typically about 30° C, and generally lower than on the dry surface). When dry and wet 
frictions are compared, the calculation shows first of all a water-induces friction decrease 
of ~ 30% at T = 60°C. The water-induced decrease becomes somewhat less (~ 20%) 
if the wet substrate temperature is (realistically) reduced to T = 30°C. We have also 
calculated /z-slip curves, and they show a similar reduction in the friction for wet road 
surfaces. These figures are in excellent agreement with the known reduction of low-speed 
rubber friction on road surfaces [53, 54]. We note in addition that the decreasing friction 
with increasing temperature shown in Fig. 42 is very commonly observed for rubber. It 
results from the shift in the viscoelastic spectrum to higher frequencies with increasing 
temperature, making the rubber more elastic and less viscous, in turn reducing the 
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rubber friction. 

The above picture does in our view catch an important novel effect of water on 
rubber friction. Yet, it is certainly open to refinements in various ways. First, dry 
friction of tires is not pure sliding but also involves some stick-slip [56] . This effect is 
included in the calculation of /x-slip curves, but the observed reduction in the effective 
friction is similar to in Fig. 42. Second, after enough time all sealing regions leak. This 
will be particularly true in the present case because the upper boundary of a water filled 
pool, which is in contact with the rubber, still has roughness on many length scales. So 
one cannot expect the rubber to make equally perfect contact everywhere, and there will 
be narrow channels through which the water slowly leaks out of the pools. As a result, 
for sufficiently low sliding velocities the negative water influence on rubber friction may 
revert to negligible. Experiments have indeed shown that for extremely low velocities 
v < 0.7 m/s the difference in ^ between dry and wet surfaces is very small [56]. We 
should also stress that the effects addressed here clearly apply only to moderately wet 
substrates and for rolling or sliding velocities v < 60km/h. For flooded surfaces and 
v > 60 km/h aquaplaning may occur, which originates instead from the inertia of the 
water. Finally, for rubber friction on relative smooth wet surfaces, where the adhesional 
interaction is important, the so called dewetting transition may be important [57, 58, 59]. 

6.4- Lubricated rubber O-ring seals 

Surface roughness has also a big influence on lubricated rubber O-ring seals [36]. Tests 
have shown that the longer a lubricated seal sits idle, the higher is static, or start-up 
friction coefficient fi s . Eventually, for smooth surfaces, the friction coefficient reaches 
a maximum almost as high as for unlubricated seal. This increase in fi s with time is 
caused by the squeeze-out of most of the lubricant from the contact area [60]. However, 
the static friction can be reduced by optimizing the surface roughness and lubricant 
viscosity. Thus, in one experiment [36] with surfaces with the rms roughness amplitude 
a PS 0.4 /zm tiny pockets with lubricant was found trapped at the interface, making it 
available at startup. Too smooth finish leaves no pockets of trapped lubricant, while 
too rough surfaces can cause high wear. 

7. Adhesion 

In this section we discuss adhesion between rough surfaces. We point out that even 
when the force to separate two solids vanishes, there may still be a finite contact area 
(at zero load) between two solids as a result of the adhesional interaction between the 
solids. We also present two applications where surface roughness influences adhesion 
in a fundamental way. We first consider an industrial application related to one step 
in the fabrication of rubber sheets. For this system a recently developed viscoelastic 
contact mechanic theory, combined with surface topography measurements using the 
Atomic Force Microscope, was able to explain the origin of a long-standing technological 
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problem. 

The second application is a study of adhesion relevant to biological systems, e.g., 
flies, crickets and lizards, where the adhesive microstructure consists of arrays of thin 
fibres and plates. The effective elastic modulus of the fibre-plate arrays can be very 
small, which is of fundamental importance for adhesion on smooth and rough substrates. 
This application illustrate how nature, through the process of natural selection, has been 
able to produce elastically very soft, but still wear- resistant layers, which can make good 
contact and exhibit strong adhesion even to surfaces with roughness on all length scales 
down to the atomic dimension. 

7.1. Adhesion between rough surfaces 

A theory of adhesion between an elastic solid and a hard randomly rough substrate 
must take into account that partial contact may occur between the solids on all length 
scales. For the case where the substrate surface is a self afnne fractal, theory shows 
that when the fractal dimension is close to 2 complete contact typically occurs in the 
macro asperity contact areas (the contact regions observed when the system is studied 
at a magnification corresponding to the roll-off wavelength Ao = 27r/g of the surface 
power spectra, see Fig. 6), while when the fractal dimension is larger than 2.5, the area 
of (apparent) contact decreases continuously when the magnification is increased. An 
important result is that even when the surface roughness is so high that no adhesion can 
be detected in a pull-off experiment, the area of real contact (when adhesion is included) 
may still be several times larger than when the adhesion is neglected. Since it is the 
area of real contact which determines the sliding friction force, the adhesion interaction 
may strongly affect the friction force even when no adhesive force can be detected in a 
pull-off experiment. 

The influence of surface roughness on the adhesion between rubber (or any other 
elastic solid) and a hard substrates has been studied in a classic paper by Fuller and 
Tabor [61] (see also [62, 63, 64, 65, 66, 67, 68]). They found that already a relative 
small surface roughness can completely remove the adhesion. In order to understand 
the experimental data they developed a very simple model based on the assumption 
of surface roughness on a single length scale. In this model the rough surface is 
modeled by asperities, all possessing the same radius of curvature, and with heights 
following a Gaussian distribution. The overall contact force was obtained by applying 
the contact theory of Johnson, Kendall and Roberts [69] to each individual asperity. 
The theory predicted that the pull-off force, expressed as a fraction of the maximum 
value, depends upon a single parameter, which may be regarded as representing the 
statistically averaged competition between the compressive forces exerted by the higher 
asperities trying to prize the surfaces apart and the adhesive forces between the lower 
asperities trying to hold the surfaces together. This picture of adhesion developed by 
Tabor and Fuller would be fine if real surfaces had roughness on a single length scale 
as assumed in their study. However, with roughness occurring on many different length 
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scales, a qualitatively new picture emerges [23], where, e.g., the adhesion force may even 
vanish (or at least be strongly reduced), if the rough surface can be described as a self 
affine fractal with fractal dimension Df > 2.5. In fact even for surfaces with roughness 
on a single length scale, the formalism used by Fuller and Tabor is only valid at "high" 
surface roughness, where the area of real contact (and the adhesion force) is very small. 
The theory that will be presented below is particularly accurate for "small" surface 
roughness, where the area of real contact equals the nominal contact area. 




Figure 43. A rubber surface is "pulled" into a cavity of the hard solid substrate by the 
rubber-substrate adhesional interaction. The elastic energy stored in the deformation 
field is of order EXh 2 . 



7.1.1. Qualitative discussion Let us estimate the energy necessary in order to deform 
a rubber block so that the rubber fills out a substrate cavity of height h and width A. 
The elastic energy stored in the deformation field in the rubber is given by 

where the stress a xs Ee, and E is the elastic modulus. The deformation field is mainly 
localized in a volume ~ A 3 (see Fig. 43) where the strain e ~ h/X. Thus we get 
U el w X 3 E(h/X) 2 = EXh 2 . 

Let us now consider the role of the rubber-substrate adhesion interaction. If the 
elastic energy U c \ ~ EXh 2 stored in the deformed rubber is smaller than the gain in 
adhesion energy U a d ~ A7A 2 , where A7 = 71 + 72 — 712 is the change of surface free 
energy (per unit area) upon contact due to the rubber-substrate interaction (which 
usually is mainly of the van der Waals type), then (even in the absence of an external 
load F N ) the rubber will deform spontaneously to fill out the substrate cavities. The 
condition U c \ = gives h/X as (A^y/EX) 1 ^ 2 . For example, for very rough surfaces with 
h/X ~ 1, and with parameters typical for rubber E = 1 MPa and A7 = 3meV/A , 
the adhesion interaction will be able to deform the rubber and completely fill out the 
cavities if A < 0.1 /im. For very smooth surfaces h/X ~ 0.01 or smaller. In that case 
the rubber will be able to follow the surface roughness profile up to the length scale 
A ~ 1 mm or longer. 
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The argument given above shows that for elastic solids with surface roughness on 
a single length scale A, the competition between adhesion and elastic deformation is 
characterized by the parameter 9 = Eh 2 /XS ~ U e \/U & d-, where h is the amplitude of the 
surface roughness and 5 = 4(1 — u 2 )A'~f/E the so called adhesion length, v being the 
Poisson ratio of the rubber. The parameter 9 is the ratio between the elastic energy 
and the surface energy stored at the interface, assuming that complete contact occurs. 
When 9 3> 1 only partial contact occurs, where the elastic solids make contact only 
close to the top of the highest asperities, while complete contact occurs when 6<cl. 

7.1.2. Pull- off force Consider a rubber ball (radius Ro) in adhesive contact with 
a perfectly smooth and hard substrate. The elastic deformation of the rubber can 
be determined by minimizing the total energy which is the sum of the (positive) 
elastic energy stored in the deformation field in the rubber ball, and the (negative) 
binding energy between the ball and the substrate at the contact interface. The energy 
minimization gives the pull-off force [69, 70] 

F c = (37r/2)i? A7. (16) 

Consider now the same problems as above, but assume that the substrate surface 
has a roughness described by the function z = h(x). Let us further assume a surface 
roughness power spectrum with a roll-off wavelength Ao = 27r/g (see Fig. 6) smaller 
than the diameter of the nominal contact area between the two solids. In this case 
we can still use the result (16), but with A7 replaced by 7 e ff. The effective interfacial 
energy 7 cff is the change in the interfacial free energy when the elastic solid is brought 
in contact with the rough substrate. 7 e fr(C) depends on the magnification (, and the 
interfacial energy which enter in the rubber ball pull-off experiment is the macroscopic 
interfacial energy, i.e., 7 e ff(C) for £ = 1. If A is the nominal contact area and A 1 the 
true atomic contact then 

A) 7cff (l) = ^iA 7 - U A (17) 

where U e \ is the elastic energy stored at the interface as a result of the elastic 
deformations necessary in order to bring the solids in atomic contact in the area A 1 . 

7.1.3. Stress probability distribution The theory in Ref. [23] is based on the contact 
mechanics formalism described in Sec. 4.1. Thus, we focus on the stress probability 
distribution function P(a,Q which satisfies Eq. (5): 

d_P = d^P 

We assume that detachment occurs when the local stress on the length scale L/( reaches 
— <7 a (C)- Thus, the following boundary condition applies to the present case 

P(-<7 a (C),C)=0 

This boundary condition replaces the condition P(0, () = valid in the absence of 
adhesion (see Sec. 4.1). 
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Let us consider the system on the characteristic length scale A = L/C,. The quantity 
<7 a (£) is the stress necessary to induce a detached area of width A. This stress can be 
obtained from the theory of cracks, where for a penny-shaped crack of diameter A 





7T7eff(C)£ 


1/2 


les(()Eq 


Oa = 






2(1 - z/ 2 )_ 



where q = 2tt/X = C<1l- In Ref. [23] we derived two equations for 7 e fr(C) and -P(C) which 
determine how these quantities depend on the magnification (; those equations are the 
basis for the numerical results presented below. 

7.1.4- Numerical results Fig. 44 shows (a) the effective interfacial energy 7 e fr(C) (C — 1) 
and (b) the normalized area of real contact, -P(Ci) = ^.(Ci)Mo, as a function of qoh . 
Results are shown for q 5 = 0.1, 0.2, 0.4 and 0.8. We will refer to 7 e ff(l) at the 
magnification ( = 1 as the macroscopic interfacial free energy which can be deduced 
from, e.g., the pull off force for a ball according to Eq. (16). Note that for qo5 = 0.4 and 
0.8 the macroscopic interfacial energy first increases with increasing amplitude ho of the 
surface roughness, and then decreases. The increase in 7 c g- arises from the increase in 
the surface area. As shown in Fig. 44(b), for small h the two solids are in complete 
contact, and, as expected, the complete contact remains to higher h as 5 ~ A^y/E 
increases. Note also that the contact area is nonzero even when 7 e fr(l) is virtually zero: 
the fact that 7 e fr(l) nearly vanishes does not imply that the contact area vanishes (even 
in the absence of an external load), but rather that the (positive) elastic energy stored 
at the interface exactly balances the (negative) adhesion energy from the area of real 
contact. The stored elastic energy at the interface is returned back when removing the 
block, and when 7 e ff(l) ~ it is just large enough to break the block- substrate bonding. 

7.1.5. Experimental manifestations of adhesion Unfortunately, the surface roughness 
power spectrum has not been measured for any of the real surfaces for which adhesion has 
been studied in detail. Instead only the roughness amplitude (center line average) and 
the radius of curvature of the largest surface asperities was determined. Nevertheless, the 
experimental data of Fuller, Tabor, Briggs, Briscoe and Roberts [61, 64, 65] are in good 
qualitative agreement with our theoretical results. In Fig. 45 we show the macroscopic 
interfacial energy for "hard" and "soft" rubber in contact with Perspex, as a function 
the substrate (Perspex) roughness amplitude as obtained by Briggs and Briscoe [64]. It 
is not possible to compare these results quantitatively with the theory developed above 
since the power spectrum C(q) was not measured for the Perspex substrate. Even if 
the surfaces would be self afline fractal as assumed above, not only would the surface 
roughness amplitude change from one surface to another, but so will the long distance 
cut off length Ao and hence also the ratio Ci — Qi/Qo- I n the experiments reported in 
Ref. [64] the Perspex surfaces where roughened by blasting with fine particles. The 
roughness could be varied through the choice of the particles and the air pressure. 

One practical problem in comparing the theory to experimental data is that most 
rubber materials have a wide distribution of relaxation times, extending to extremely 
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Figure 44. (a) Macroscopic interfacial energy for adhesion of rubber to a fractal 
surface with H = 0.8 and qi/qo = Ci = IOO7 as a function of qoh . (b) Normalized 
area of real contact, P((i) = A((i)/Ao, as a function of qoho. For qoS = 0.1, 0.2, 0.4 
and 0.8 as indicated. 



long times. This effect is well known in the context of rubber friction (see Sec. 6.1), where 
measurements of the complex elastic modulus show an extremely wide distribution of 
relaxation times, resulting in large sliding friction even at very low sliding velocities, 
v < 10~ 8 m/s. 

The effect of the stored elastic energy on adhesion has recently been studied using a 
polyvinylsiloxane rubber block squeezed against a smooth glass surface for a fixed time 




Roughness amplitude (urn) 

Figure 45. The macroscopic interfacial energy (obtained from the pull-off force) for 
a smooth rubber surface (ball) in contact with a Perspex surface as a function of the 
roughness (center line average) of the Perspex. Results are shown for a "soft" rubber 
(E = 0.063 MPa) and a "hard" rubber (E = 0.487 MPa). From [64]. 

period before measuring the pull-off force [71] . The square-symbols in Fig. 46 show the 
pull-off force as a function of the squeezing force. For squeezing forces F N > 850 mN the 
pull off force decreases. This may be explained by a drastic increase of the elastic energy 
stored in the rubber because of the strong deformation of the rubber, see Fig. 47(top), 
some of which remains even when the load is removed as a result of the rubber-glass 
friction at the interface. This energy, freed during the process of unloading, will help to 
break the adhesive bonds at the interface. This effect is even stronger when the surface 
is structured. Thus, the triangles in the figure shows the pull-off force when the rubber 
surface is covered by a regular array of rubber cylindrical asperities. In this case the 
pull-off force drops to nearly zero for F N > 700 mN. Visual inspection shows that in this 
case the cylindrical asperities at high load bend and make contact with the glass on one 
side of the cylinder surface, see Fig. 47(bottom). This again stores a large amount of 
elastic energy at the interface which is given back during pull-off, reducing the pull-off 
force to nearly zero. 

7.1.6. The role of plastic yield on adhesion When the local stress in the asperity 
contact regions between two solids becomes high enough, at least one of the solids yields 
plastically. This will tend to increase the effective adhesion (or pull-off force) for the 
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Figure 46. The pull-off force as a function of the squeeze force or load, for silicon 
rubber in contact with a smooth glass surface. From Ref. [71]. 



following three reasons. First, the area of real contact between the solids will increase as 
compared to the case where the deformations are purely elastic. Secondly, the amount 
of stored elastic energy in the contact regions (to be given back at pull-off) will be 
reduced because of the lowered elastic deformations. Finally, for many materials plastic 
yield will strengthen the junctions [72]. For example, most metals are protected by thin 
oxide layers, and as long as these are intact the main interaction between the surfaces 
in the contact areas may be of the van der Waals and electrostatic origin. However, 
when plastic yield occurs it may break up the oxide films resulting in direct metal-metal 
contact and the formation of "cold- welded" junctions. When this occur, because of the 
high ductility of many metals, during pull-off "long" metallic bridges may be formed 
between the solids so that instead of having junctions popping one after another during 
pull-off, a large number of adhesive junctions may simultaneously impede the surface 
separation during pull-off, leading to a large pull-off force. However, experiment have 
shown [73] that just squeezing before pull-off will in general only result in very few cold 
welded junctions, while squeezing and sliding will break up the oxide film, resulting in 
the formation of many more cold welded contact regions, and will hence result in a much 
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Figure 47. Elastic deformation of a rubber block with a smooth surface (top) and 
a structured surface (bottom), (a) shows the initial state before applying a squeezing 
force, and (b) the new state (without load) after applying (and then removing) a very 
large squeezing force. In state (b) a large elastic energy is stored in the rubber which 
is "given back" during pull-off resulting in a nearly vanishing pull-off force. 

larger pull-off force. 

7.2. The adhesion paradox 

The biggest "mystery" related to adhesion is not why it is sometimes observed but 
rather why it is usually not observed. Even the weakest force in Nature of relevance 
in condensed matters physics, namely the van der Waals force, is relatively strong on 
a macroscopic scale. For example, even a contact area of order 1 cm 2 could sustain the 
weight of a car (i.e., a force of order 10 4 N) [see Fig. 48(a)] even if only the van der 
Waals interaction operated at the interface. [Here we assumed that the bond breaking 
occur uniformly over the contact area as illustrated in Fig. 48(b).] However, this is never 
observed in practice and this fact is referred to as the adhesion paradox. 

There are several reasons why adhesion is usually not observed between macroscopic 
bodies. For example, on a macroscopic scale the bond-breaking usually does not occur 
uniformly as in Fig. 48(b), but occurs by crack propagation, see Fig. 48(c). The local 
stress at the crack tip is much higher than the average stress acting in the contact area, 
and this drastically reduces the pull-off force. Another reason, already addressed in 
Sec. 7.1, is the influence of surface roughness. Thus, for elastically hard surfaces the 
true atomic contact between the solids at the interface is usually much smaller than the 
nominal contact area. In addition, the elastic energy stored in the solids in the vicinity 
of the contact regions is given back during pull-off and helps to break the interfacial 
bonds between the solids (see Sec. 7.1). 

It is interesting to note that for very small solid objects, typically of order 100 fxm 
or smaller, the bond breaking may occur uniformly over the contact area (no crack 
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propagation) so that adhesion between smooth surfaces of small objects, e.g., in 
micromechanical applications (MEMS), may be much stronger than for macroscopic 
bodies, and this fact must be taken into account when designing MEMS [74, 75]. 
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Figure 48. Even the weakest force in Nature which is of relevance in condensed 
matters physics, namely the van der Waals force, is relative strong on a macroscopic 
scale. Thus, for example, if the bond breaking were to occur uniformly over the 
contact area as in (b), already a contact area of order 1 cm 2 could sustain the weight 
of a car (i.e., a force of order 10 4 N) [see (a)]. However, on a macroscopic scale the 
bond-breaking does not usually occur uniformly over the contact area, but by crack 
propagation, see (c), which drastically reduce the pull-off force. Secondly, interfacial 
surface roughness drastically reduces the pull-off force. Finally, the stored elastic 
energy will largely balance the adhesion energy (see text) 



7.3. Adhesion in rubber technology 

Most rubber compounds of commercial use contains a large fraction (around 30%) of 
filler particles which usually are mixtures of carbon and silica particles. The silica 
particles are very small (below micrometer-size) and very hard. As a result, when 
rubber is sliding on a substrate, even if the latter is very hard, e.g., stone or steel, the 
substrate will get polished by the rubber. We have discussed this effect above for road 
surfaces (see Sec. 6.2). Here we describe another recent observation of the same effect 
with important practical implications. 

One stage in the production of rubber for technological applications involves the 
mixing in open mills, where unvulcanized rubber (with filler particles) is fed between 
rotating steel cylinders. The slip of the (silica-particle containing) rubber relative to 
the steel walls during the mixing process result in a continuous "polishing" (or wear) 
of the steel surfaces. This will slowly increase the rubber-steel contact area, and finally 
the rubber may adhere to the steel surfaces, which is of course unwanted. 
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Using the Atomic Force Microscopy (AFM) we measured the surface topography of 
the steel cylinders, and found that during production the steel surfaces are continuously 
being polished by the rubber containing silica particles. This leads to a continuously 
increased contact area between the steel and the rubber sheet, and to an increased 
adhesive steel-rubber interaction. Using a recently developed viscoelastic contact 
mechanics theory (see Sec. 4.3) we further found that the increase in the contact area 
can be very large. 

Also shown by AFM measurements was that lapping of the steel surface by 
corundum paper increases the steel surface roughness (see Fig. 49), which is easy to 
understand when one recognizes that the corundum particles have on the average larger 
diameter than the silica particles in the rubber. Thus we conclude that lapping or 
"polishing" may increase the roughness of the steel cylinders, and hence reduce the 
adhesive interaction to such a level that it becomes unimportant. 
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Figure 49. The surface roughness power spectra of a steel cylinder before and after 
lapping with corundum paper. The lapping increases the surface roughness. 



7.4- Adhesion in biology 

How can a fly or a cricket walk on a glass window, or a lizard (see Fig. 50) move on a stone 
or concrete wall? These fundamental questions have interested children and scientists 
for many years, and recently very important experimental work has been performed 
that gave a deeper insight into these questions [76, 77]. Here we focus mainly on dry 
adhesion, which seems to be relevant for lizards, and we discuss the influence of surface 
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Figure 50. A lizard adhering to a nearly vertical surface. 

roughness on the adhesion between a lizard or a gecko toe and a rough hard substrate. 

It has been demonstrated that a foot of a gecko can adhere to a substrate with 
a force ~ ION (corresponding to the weight of 1 kg!). The typical weight of a Tokay 
gecko lizard (Gecko gecko) is approximately 40 g meaning that only 1% of the maximum 
adhering force of its feet is required to support the whole weight of the gecko. This raises 
the question of why geckos are apparently so over-built. However, a gecko must be able 
to adhere to very rough surfaces, and we will show below that the adhesion to rough 
surfaces can be reduced significantly [78, 79]. 

The adhesion between two bodies in contact results almost entirely from the regions 
where the surfaces are separated by one nanometer or less. For hard solids this area of 
real contact is extremely small. Furthermore, for elastically hard solids a large elastic 
energy is stored in the solids in the vicinity of the contact regions and during separation 
of the solids this energy is released and will hence reduce the separation force as explained 
earlier. Thus, strong adhesion is only possible if at least one of the solids is elastically 
very soft, or if there is an elastically very soft layer at the interface between the solids. 
Since the lizard skin comprises a relative stiff material (keratin), with an elastic modulus 
of the order 1 GPa, i.e., about 1000 times higher than tire rubber, it is not immediately 
obvious why the lizard can adhere to very rough stone walls. 

As it turns out, during millions of years of evolution and natural selection, an 
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extremely soft elastic layer has appeared on the lizard pad surface. This layer is built in a 
hierarchical manner from fibres and plates (Figs. 51 and 52), that reflect the hierarchical 
nature of most natural surfaces (to which the lizard must be able to adhere), which have 
roughness on all length scales, from the macroscopic scale (e.g., the size of the lizard 
toe pad) down to the atomic scale. The skin of the lizard pad, which consist of a 
~ 100 fim thick keratin layer, is covered by a dense layer of fibres or hair (setae) (length 
rs 100 yum and thickness ~ 4/iin). Each of these fibres branches out into about 1000 
thinner fibres (length ~ 10 jum and width ~ 0.1 /im), and each terminal fibre ends with 
a thin (5 — 10 nm) leaf-like plate (spatula). This construction makes the lizard adhesive 
system elastically very soft on all relevant length scales (from mm to nm). 



skin on foot pad 




Figure 51. Schematic picture of the lizard adhesive system. The skin of the lizard is 
covered by a dense layer of thin fibres or hair (setae) (length w 100 /im and thickness 
of fibre of order <~ 4/im). Each of these fibres branches out into about 1000 thinner 
fibres (length ~ 10 /im and width of order ~ 0.1 /im). Each of the thin fibres ends with 
a thin (5 — 10 nm) leaf- like plate (spatula). 

The skin of the gecko toe-pad is able to deform and follow the substrate roughness 
profile on length scales much longer than the thickness d ~ 100 /im of the elastic keratin 
film, say beyond ~ 1000 /mi. At shorter length scales the keratin film, because of its high 
elastic modulus (of order 1 GPa), can be considered rigid and flat. Elastic deformation 
of the pad surface on length scales shorter than ~ 1000 /im, involves the compliant setae 
fibre array system, with fibres of thickness ~ 4 /jm. In Ref. [78] one of us has shown 
that if the surface roughness root-mean-square amplitude, measured over a patch D x D 
with D xs 1000 /im, is smaller than a characteristic length (the adhesion length) (see 




Figure 52. Details of attachment system of the Tokay gecko (Gecko gecko). A. 
Scanning electron microscopy (SEM) micrograph of setae (st) located on thin keratin 
film. B. Magnification (SEM micrograph) of the area surrounded by the white 
rectangle in A, showing terminal branches (tb) of setae with the spatula (sp). C. 
Transmission electron microscopy micrograph of ultrathin section of two terminal 
branches (tb) with spatulae (sp). From [79]. 

Ref. [78]), then the fibre array system is able to deform (without storing much elastic 
energy) to follow the surface roughness in the wavelength range 10 < A < 1000 /xm. 
However, if the setae fibre tips were instead blunt and compact, they would not be 
able to penetrate into the surface "cavities" with diameters less than a few pim. Thus, 
negligible atomic contact would occur between the surfaces, and adhesion would be 
negligible. For this reason, there is an array of ~ 1000 thinner fibres (diameter of order 
~ 0.1 /im) at the tip of each long (thick) fibre. These fibres are able to penetrate into 
the surface roughness cavities down to length scales of a few tenths of a micrometer, see 
Fig. 54. Again if the thin fibres had blunt and compact tips made of the same "hard" 
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Figure 53. Spider man. Somewhere in our Universe there may be a planet with 
'animals' as large as humans moving on vertical walls. . . 'SPIDER-MAN'. Motion 
Picture © 2002 Columbia Pictures Industries, Inc. Spider-Man Character ® & © 
2002 Marvel Characters, Inc. All Rights Reserved. Courtesy of Columbia Pictures. 

keratin as the rest of the fibre, one would still obtain very small adhesion, since much 
elastic energy would be needed to deform the surfaces of the thin fibres to make atomic 
contact with the substrate. But in fact the top of the thin fibres end with thin leaf-like 
plates, which can easily bend (without storing much elastic energy) to follow the surface 
roughness profile [79, 80]. The calculations presented in Ref. [79] show that, for rough 
surfaces with the fractal dimension D{ > 2.3, very small spatula-substrate adhesion 
may occur in most cases. However, natural surfaces tend to have a fractal dimension of 
order 2.2, and adhesion may be appreciable even for very rough surfaces in these cases. 
Experiments to test the theoretical results are in progress. 

There is an interesting aspect of fibre adhesion which differs from that involving 
compact solids. Fibre adhesion depends only on the height probability distribution Ph 
and not on the power spectra C(q). The reason is that the elastic energy stored in a 
bent fibre depends only on the separation between the surfaces at the interface which 
is determined by the surface roughness height distribution P h . Thus, the elastic energy 
stored in a bent fibre does not depend directly on C(q). On the other hand for compact 
solids the elastic energy is determined only by C(q). 

Lizards are the heaviest living objects on this planet that are able to adhere to, 
for example, a rough vertical stone wall. Since the surface area of a body increases less 
than the volume, or mass as the size of the body increases, the adhesive system in large 
living bodies, such as lizards, must be much more effective (per unit attachment area) 
than in smaller living objects such as flies or beetles. This implies that lizards have 
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Figure 54. The Lizard fibre-arid-plate adhesive system can be easily elastically 
deformed to bind to a rough substrate, even when the substrate has roughness on 
all length scales from, say, ~ 1 mm to ~ 1 nm. 

the most effective adhesive systems found in the biological evolution for the purpose of 
locomotion. This is confirmed by electron microscopy studies. Thus, the spatula are 
thinner in lizards than in beetles. Also the diameter of the terminal branches is smaller. 
This implies that less elastic energy per unit surface area will be stored in the lizard 
adhesive system, resulting in a stronger adhesion for lizards than for beetles. 

The construction of man-made adhesives based on fibre and plate arrays might 
constitute an attractive alternative to the usual adhesives based on thin polymer films 
(see Fig. 53 for one "application"). Some pioneering experiments have indeed shown 
enhanced adhesion for fibre array systems, but no man-made systems with the hierarchic 
nature found in biological systems have so far been produced. In addition, if the fibres 
are too long and thin, or the fibre material is (elastically) too soft, the attractive 
fibre-fibre interaction will result in elastic instabilities, leading to fibre bundles or fibre 
"condensation" into compact layers. This effect was in fact observed in the first man- 
made fibre adhesion system (which used fibres made from rubber which is 1000 elastically 
softer than keratin used in most biological applications), see Fig. 55, and was at the 
same time predicted theoretically [78]. It turns out that the fibre arrays in Lizards 
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are close to this instability, but, unsurprisingly, on the correct side of it, and no fibre 
bundling or condensation occur in these biological systems. 




Figure 55. When the fibres are too thin and long, or made of materials clastically too 
soft (as in the present case of rubber shown here) the standing-up fibre array system in 
(a) is unstable, and after contact with a substrate the fibres end up in the bundled-up 
state (b). Adapted from [81]. Reproduced with permission. 

7. 5. The role of liquids in adhesion between rough solid surfaces 

As explained in Sec. 7.1, surface roughness reduces the adhesion between clean surfaces. 
First, it lowers the area of real contact. Since the adhesion interaction comes almost 
entirely from the area where the solids make atomic contact, it is clear that the surface 
roughness can drastically reduce the adhesion. Secondly, elastic deformation energy is 
stored in the vicinity of the asperity contact regions. During pull-off the elastic energy 
is "returned back" to the system, usually resulting in a drastic reduction in the effective 
adhesion and the pull-off force. 

Most surfaces have at least nano-scale roughness, and hard solids in the normal 
atmosphere have at least a monolayer of liquid-like "contamination" molecules, e.g., 
water and hydrocarbons. Small amounts of a wetting lubricant or contamination liquids 
between rough solid walls may drastically enhance the adhesion. Thus, for surfaces with 
nanoscale roughness, a monolayer of wetting liquid may result in the formation of a large 
number of nano-bridges between the solids, which increases the pull-off force. This effect 
is well known experimentally. For example, the adhesion force which can be detected 
between gauge blocks (steel blocks with very smooth surfaces) is due to the formation of 
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many very small capillary bridges made of water or organic contamination. For thicker 
lubrication or contamination films the effective adhesion will be more long-ranged but 
the pull-off force may be smaller. The thickness of the lubricant or contamination layer 
for which the pull-off force is maximal will in general depend on the nature of the surface 
roughness, but is likely to be of order the root-mean-square roughness amplitude. In 
fact, it is an interesting and important problem to find out at exactly what liquid 
thickness the pull-off force is maximal. 

Some insects such as flies or crickets inject a thin layer of a wetting liquid in the 
contact region between the insect attachment surfaces and the (rough) substrate. The 
optimum amount of injected liquid will depend on the nature of the substrate roughness, 
and it is likely that the insect can regulate the amount of injected liquid by a feedback 
system involving the insect nerve system. 

Here we consider the adhesion between two solid elastic walls with nanoscale 
roughness, lubricated by octane [68, 60, 67]. We consider two types of substrates (bottom 
surface) - flat and nano-corrugated (corrugation amplitude 1 nm and wavelength of the 
corrugation in x and y direction, 4nm) - and varied the lubricant coverage from ~ 1/8 
to ~ 4 monolayers of octane. The upper surface (the block) is assumed to be atomically 
smooth but with a uniform cylindrical curvature with a radius of curvature R ~ 100 nm 
(see Fig. 58 below). The simulation results presented here were obtained using standard 
molecular dynamics calculations [67]. 

Fig. 56 shows the variation of the average pressure during retraction as the block 
moves a distance of 16 A away from the substrate. The pull-off (retraction) velocity 
was v z = 1 m/s. We varied the lubricant coverage from to 1 monolayer in the contact 
region. The pull-off force is maximal when the adsorbate coverage is of the order of one 
monolayer [curve (f)]. However, the pull-off force is still smaller than for a flat substrate 
without lubricant [curve (a)]. As a function of the octane lubricant coverage and for the 
corrugated substrate, the pull-off force first increases as the coverage increases from zero 
to ~ 1 monolayer, and then decreases as the coverage is increased beyond monolayer 
coverage (not shown). 

At low octane coverage, the octane molecules located in the substrate corrugation 
wells during squeezing are pulled out of the wells during pull-off, forming a network 
of nano capillary bridges around the substrate nanoasperities, thus increasing adhesion 
between two surfaces, see Figs. 57 and 58. For greater lubricant coverages a single 
capillary bridge is formed. 

Let us discuss the nature of the adhesion for the corrugated substrate, with about 
1/4 monolayer of octane in the contact region. Fig. 57 shows snapshot pictures of the 
lubricant layer during retraction, as the block moves away from the substrate for three 
different block positions d — 0, 3 and 6 A. Only the central part of the contact between 
the block and the substrate is shown, top view, after removing the block and substrate 
atoms. In the beginning (d = OA) octane molecules are located in the substrate 
corrugation wells, or cavities with direct metal-metal contact between the block and 
the top of the substrate nano asperities (see Fig. 58). During retraction (d — 3 A) the 
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distance (Angstrom) 

Figure 56. Simulated adhesion between two solid elastic walls, one of which with 
nanoscale roughness, and with a variable amount of wetting lubricant in between. 
The variation of the average pressure during retraction develops as a block moves a 
distance of 16 A away from the substrate. Octane CsHig was used as lubricant. Pull- 
off (retraction) velocity was v z = lm/s. (a) For the flat substrate without lubricant, 
(b) For the corrugated substrate without lubricant. Curves (c)-(f) show results for 
the corrugated substrate with about 1/8, 1/4, 1/2 and 1 monolayer of octane in the 
contact region, respectively. For clarity, the curve for the flat substrate (a) is displaced 
to the right, by 2 A. 

octane molecules are pulled out of the wells forming an almost symmetric network of 
nano-bridges around the asperity tops, increasing the adhesion between the two surfaces. 
This configuration corresponds to the maximal adhesion force, see curve (d) in Fig. 56. 
Thus maximal adhesion is achieved via the formation of many small capillary nano- 
bridges, involving just a few molecules for each bridge (see Fig. 58). Further retraction 
(d — 6 A) results in the collapse of the nano-bridges and the formation of a single "large" 
capillary bridge in the centre of the contact region. 

8. Summary and outlook 

Surface roughness has a huge influence on many common phenomena. It is the main 
reason for why macroscopic bodies usually do not adhere to each other with any 
measurable strength. For example, if the floor and the sole of the shoes were atomically 
smooth and clean, it would be impossible to walk. The (near) absence of adhesion 
in most situations is crucial for the function of many man-made and nature-made 
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d = 6 A 




Figure 57. Snapshot pictures (for three different block positions d = 0, 3 and 6 A) 
of the lubricant layer during retraction. We only show the lubricant molecules in the 
central part of the contact area between the block and the substrate surfaces (top 
view, surfaces parallel to the plane of the image). Corrugated substrate with about 
1/4 monolayer of octane in the contact region. The circles indicate the position of 
several asperity tops of the corrugated substrate surface. 
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d = OA 




d = 10 A 




Figure 58. Snapshot pictures (for six different block positions) during retraction. The 
snapshot pictures show the side view of the central 108 A x 50 A section (in the icy- 
plane) of the contact area. Octane CsHis was used as lubricant. Pull-off (retraction) 
velocity was v z = lm/s. For the corrugated substrate with about 1/4 monolayer of 
octane in the contact region. 

constructions. 

The surface to volume ratio of solid objects increases as the lateral size of the object 
decreases. The role of surface roughness becomes therefore more important as the size of 
objects decreases. The present drive toward the miniaturization of mechanical devices, 
e.g., Micro-Electro- Mechanical Systems (MEMS), requires a better understanding of the 
role of surface roughness on, e.g., contact mechanics and adhesion. 

Surface roughness is also of great importance for the function of many biological 
systems. Flies, bugs, crickets and lizards have developed very soft layers on their feet 
organs which allow them to attach and move on both very smooth and rough vertical 
solid walls, e.g. stone walls or leafs. Another example is non-wetting coatings on plant 
surfaces based on surface roughness on many different length scales (the so called Lotus 
effect [3]). 
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The roughness of surfaces can nowadays be studied rather straightforwardly using 
standard equipments based on optical methods and on cantilever methods, e.g., the 
atomic force microscope (AFM). These methods cover the whole length scale from 
atomic dimension to macroscopic distances. Thus, the AFM can probe the surface profile 
from ~ 1 nm to 100 /iin and optical methods from ~ 1 fim to kilometers. For randomly 
rough surfaces, the most important quantity which can be deduced from the measured 
height profile is the surface roughness power spectrum. We have shown throughout 
this paper how the roughness power spectrum determines the contact mechanics and 
adhesion for solid objects in direct contact. It also governs rubber friction on rough 
substrates, e.g., friction of tires on a road surfaces, and influence other phenomena of 
technological importance, e.g., the roughness induced leaking of sealings. 

Thus study and characterization of surface roughness is important not only for 
understanding many natural and biological phenomena, but also for many technological 
processes. The present drive toward miniaturization and the design of optimal systems 
by transfer of ideas from biology to materials science (bionics) [3, 82], is likely to 
accelerate the interest and efforts to study and predict the influence of surface roughness 
on many phenomena. 
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Appendix A. More about surface roughness 

Assume that a randomly rough surface is described by the height profile z = h(x), 
where x = (x, y) is a 2D vector in the surface plane z = 0, conveniently chosen so that 
(h) = and so that ((h — (h)) 2 } is minimal. The statistical properties of randomly 
rough surfaces is completely specified when all the correlation functions 

(h(x 1 )h(x 2 )), (h(x 1 )h(x 2 )h(x 3 )) , 

are known. Here (. . .) stands for ensemble averaging. In most cases one assume that 
the correlation functions involving an odd number of ^-functions vanishes, while the 
correlation functions involving an even number of /i-functions can be decomposed into 
a product of pair correlation functions, e.g., if we denote h(xi) = hi and so on, 

{hih 2 hzh A ) = {hih 2 ){h 3 h A ) + {hih z ){h 2 h A ) + {h-Ji^) (h 2 h 3 ) 

In this case the surface is completely specified by the surface roughness power spectrum 
C(q), which is the Fourier transform of the pair correlation function (h(x)h(o)), and 
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the surface height distribution P h is Gaussian (see below). 

It is possible to generate surface roughness profiles, which are very similar to 
experimentally observed surface profiles, as follows: The surface height over alxl 
square area can be expressed through its Fourier series: 

h( x ) = Y J B( q y [q - x+m] (a.i) 

where q spans all the vectors whose components are whole multiples of 2tt/L. Since h(x) 
is real, B(—q) = B(q) and <f>(—q) = —0(g). If 4>{q) are independent random variables, 
uniformly distributed in the interval [0, 2ir[, then one can easily show that higher order 
correlation functions can be decomposed into a product of pair correlations in the way 
described above. Here we will demonstrate that in this case the height probability 
distribution Ph is always Gaussian, while any power spectrum can be arbitrarily imposed 
by choosing properly the amplitudes B(q) 

Let us consider a randomly rough surface described by equation (A.I). Provided 
that the phases <p(q) are uniformly distributed and independent, the statistical 
properties of the surface are translationally invariant, thereafter 

(h(x 1 )h(x 2 )) = C(xi - x 2 ). 
The surface roughness power spectrum is defined by 

C (Q) = 7^ / d " x C(x)e- i<iX (A.2) 



(2tt) 2 „ 

By substituting (A.I) in (A.2) and using 
(e W) e W)> = 5q , 



-q" 



it follows 



C(q) = p/d 2 ,E |3(q')|¥^>* = £ \B{q')M<l ~ <?') 



q' q' 

When the sampling of the qr-space is dense enough we can approach the continuous limit 
by replacing 

where A is the nominal surface area. This gives 
CM = j^\B(^ 

Thus, if we choose 

B(q) = (2 7 r/L)[C(q)} 1 / 2 , 

where L = A 1 / 2 , then the surface roughness profile (A.I) has the surface roughness 
power density C(q). We can also guarantee that the statistical properties of the rough 
surface are isotropic by imposing B(q) = B(q), then C(q) = C(q) is a function of 
the magnitude q — \q\, but not of the direction of q. (The condition above shows the 
proportionality relation between the Fourier transform of the height correlation and the 
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amplitude of the corresponding Fourier component of the surface profile. Such result, 
known as Wiener-Khintchine theorem, is further discussed in Appendix C) 

Let us prove that surfaces whose correlation functions (hih 2 • • • h n ) vanish for odd 
n while for even n can be decomposed into a product of pair correlation functions, 
have Gaussian height probability distributions. First note that the height probability 
distribution 

P h = (5[h - h(x)}) = -L / da (e 1 ^-^) 

da e iah (e- ia *^) (A.3) 



1 

2^ 



But 



nl 



However, ([h(x)] n ) vanishes for odd n. Thus, writing n = 2m (m = 0, 1, ...) gives 

(—in) 2171 

<e-^ W > = E i (2^r<[M*)] am > (A-4) 
We will now make use of the following equation 

([M*)] 2m > = { ^( h2 r ( A - 5 ) 

The prefactor (2m)!/m!2 m is easy to understand: when decomposing (hih 2 ■ ■ . h 2m ) 
there are (2m)! possible ordering of hi, h 2 , ■ ■ ■ , h 2m . However, many of these products 
of pair correlation functions are identical. Thus, since in the pair correlation function 
(hihj) the order of hi and hj is irrelevant, we must divide (2m)! with 2 m . Furthermore, 
permutation of the m different pair correlation functions in a product term results in 
m! identical terms, which is the origin of the 1/m! factor in (A. 5). Substituting (A. 5) 
into (A. 4) gives 

(e -iaM*)\ = y tl^ln h 2 )r = e -§>^> 

Substituting this result into (A.3) and performing the integration over a gives 

p _ 1 c -feW 
h (2vr)VV 

where a is the rms roughness amplitude, a 2 = ([h(x)] 2 ). 

Fig. Al shows a randomly Gaussian random surface generated with the method 
discussed above, through formula (A.l). The power spectrum is imposed to be 
proportional to g~ 3 - 6 for large wave vectors, so that the Hurst exponent is 0.8 and 
the fractal dimension is 2.2. 



Appendix B. Hurst exponent and fractal dimension 



A surface is said to be self-similar if it has the same statistical properties as a magnified 
version of itself. Thus, if z = h(x,y) is the equation of a self-similar surface, then its 
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Figure Al. Computer generated rough substrate. The surface has root mean square 
roughness 0.001m, roll-off wave vector q n — 1000m -1 , and Hurst exponent H = 0.8. 
Length in the picture are in meters. 



magnified version z = \h(x/\,y/\) cannot be distinguished from the original. A self- 
affine surface is analogous, except that the magnification factor along the z direction 
differs from the in-plane magnification factor A. In particular the Hurst exponent H 
defines the scaling factor \ H along the z direction in order to recover the same statistical 
properties. The transformed surface z = \ H h(x/\,y/\) "looks" exactly as the original 
one z = h(x,y). The Hurst exponent is in between and 1, the latter being the case of 
a self-similar surface. 

The general form of the power spectrum for an isotropic self-afnne fractal surface 
is easy to derive. We have 

C(q) = J^y 2 jd 2 x (h(x)h(o))e-^ x 

Now let us write x = x' / A so that 

C(g) = j^y 2 J dV \~ 2 (h(x'/\)h(o))e-^'^ 

But since for a self affine fractal surface 

(X H h(x'/X)X H h(o)} = (h(x')h(o)) 

we get 

C(q) = T^ 2 J dV \- 2 ~ 2H (h(x')h(o)}e-^'^ 
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Thus, if we choose A = q and denote q = q/q we get 

C (q) = g-aa+H) JL / d V (h(x')h(o))e-^' 

Thus the power spectrum of a self-affine surface decreases as g" 2 ( H+1 ) with increasing 
wave vector q = \q\. 

The Hurst exponent is directly related to the fractal dimension Df through the 
formula Df = 3 — H . The proof is straightforward: consider a patch of the surface with 
extension L x L in the two in-plane directions x and y. The fractal dimension can be 
defined through the number of cubes of size A required to cover completely the surface: 
N(A) oc A~ Df for A — > 0. To cover the LxL area a number of cubes (L/A) 2 is required, 
that is the fractal dimension cannot be smaller than 2. In any A x A sub-domain the 
^-coordinate of the surface spans a range of values Az proportional to A H , because of 
the self-affine property. Actually Az = K(A/L) H , where K is the range of z values 
spanned by h(x, y) over the whole LxL domain. Provided that A is small enough, the 
corresponding Az gets larger than A itself, since H < 1. Thereafter in any sub-domain 
of size A we have to employ K(A/L) H /A oc A^ -1 cubes, and the total number of cubes 
necessary to cover the surface is thus proportional to A~( 3 ~ H \ i.e. 

Df = 3-H. 



Appendix C. Moments of power spectra 

We consider rough surfaces for which the statistical properties are isotropic and 
translational invariant. Thus the surface roughness power spectrum 

C(q) = J^y 2 J d2x (M* + x')h(x'))e-^ x (C.l) 

is independent of x' and of the orientation of the wavevector q. If we express h(x) 
through its Fourier transform 

h(x) = J d 2 q h{q)e iq - x (C.2) 

and we substitute it into (C.l) we get 

C(q) = j^y 2 J d 2 q'J d 2 q" J d 2 x (^(qO^lg"))^^^^^^^")-' 

Generally h(q) is not an ordinary function but is a distribution, thus (h(q')h(q")) 
can diverge. However we can replace h(q) with ^a(q), the Fourier transform of h(x) 
restricted to a square area A, i.e., 

1 r+y/A/2 r+VA/2 

hA{q) = j—r^ I dxi / dx 2 h(x)e iqx , 

(2tt) 2 J-VA/2 J-VA/2 

and the equations above still hold in the limit A — > oo. 



C{q) = J^y 2 j d 2 q'J d 2 q" J d 2 x (h A (q')h A (q"))e 
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To simplify this equation we use the standard relation 

J d 2 x j q x = (2n) 2 5(q) (C.3) 

to get 

C(q) = J dY {h A (q)h A (q"))jW>* (C.4) 

This expression is independent of x', at least in the limit A — > oo, hence (h(q)h(q")) = 
lim A ^, 00 {h A (q)h A (q")) must vanish when q + q" ^ o. Moreover the isotropy implies 
that (h(q)h(—q)) depends on q only through its modulus q; thereafter: 

(h(q)h(q')) = C(q)5(q + q') (C.5) 

Since (C.4) is independent of x' we may integrate over x' and divide by the surface area 
A to get 

C(q) = I / dV/ d 2 q" (h A (q)h A (q"))jW> x ' 



(27T) 



2 



(h(q)h(-q)), 



A 

where we have again used (C.3). Since h(x) is real, h A (—q) = h* A (q), and we can write 

C(q) = ^(\h A (q)\ 2 ) (C.6) 

Equations (C.5) and (C.6) are very important and useful equations of general validity. 
As an application, let us calculate the average surface slope. Using (C.2) we get 

((V/i) 2 > = J d 2 qd 2 q' (-q ■ q'){h(q)h(q')}e^ +(1 '> x 
Using (C.5) this gives 

((V/i) 2 ) = J d 2 q q 2 C{q) = 2it J dq q 3 C{q) 
For a self affine fractal surface C(q) ~ q~ 2 ( 1+H ) and 

/(T7U\2\ f 1 A 1-2H 2(1-//) 2(1- H) ^ 2(1~H) 

"'go 

if gi ^> go- Thus the average slope (and the average curvature) is determined by the 
shortest wavelength roughness components, as indicated in Fig. CI. On the other hand, 
as shown in Sec. 2, the rms-roughness amplitude is determined by mainly by the longest 
surface roughness wavelength components, i.e., the region around q ~ q (see Fig. CI). 

Appendix D. Numerical recipes for calculating power spectra 

Here we describe how the roughness power spectrum C(q) can be numerically determined 
from the surface height h(x) measured over a square area < x < L and < y < L. 
The surface roughness power spectrum is given by 



C(q) = J^yJd 2 x (h( X )h(o))e-^, 
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Figure CI. Surface roughness power spectra of a surface which is self affine fractal 
for q 1 > q > q . The rms roughness amplitude, and the average slope (and the average 
curvature), are determined mainly by the encircled regions of the power spectrum. 



which also can be written as (see Appendix C) 

C(q)= ( ^(\h A (q)\ 2 }, 
where A = L 2 is the surface area under study and 
h A (q) = j±y 2 j^x h(x)e-^. 

In order to calculate C(q) numerically from the measured height profile h(x), we assume 
that the measurement points form a square mesh with the lattice constant a, and with 
N data points along the x and y-coordinates. In other words, h(x) is sampled and its 
values are known only in the points 

x — ( n x, n y )a = na = x n , 

where n x = 1, 2, . . . , N and n y = 1,2, ... ,N are integers. The integral for the Fourier 
transform h\ can be approximated by a discrete sum: 

M<z) » -f^J2 h ^~ i{qxnxa+qynya) 

where h n = h(x n ). The Fourier transform can be sampled too, by using a grid with 
lattice spacing qi = 2tt/L = 2ir/(Na), since the finite size of the surface is limiting the 
resolution in the q-space. Thus q = (q x ,q y ) = (2nm x /L, 2irm y /L) where m x and m y 
are integer numbers between and N — 1 (any larger value of m x or m y corresponds to 
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wavelengths shorter than a, and does not carry extra information because of the finite 
resolution in the real space). Now we can write: 



2 2 

a , _;27r/™ „ i jyvi „ ^ a 



Mflf) « Tf^E^e-^^- = j^Hn, (D.l) 

where H m is the 2-dimensional Discrete Fourier Transform of h n . The big advantage 
of evaluating the power spectrum through the Fourier transform of h(x), instead of 
carrying on the calculation of the correlation, is that there exists an efficient algorithm 
for the Discrete Fourier Transform: the well known Fast Fourier Transform (FFT) 
method. Indeed the computational time scales as A 2 log(A) (in two dimensions), 
while the convolution integral to estimate the 2-dimensional correlation would imply 
a computational time that scales as iV 4 . 

A crude result can be obtained by substituting (D.l) into (C.6): the power spectrum 
obtained so far would be identical to the one resulting from the Fourier transform of 
the correlation, provided that the correlation is calculated applying periodic boundary 
conditions over the surface of area A. But real surfaces are portions of larger samples, 
where there is no match at the boundary: the periodic boundary conditions introduce 
discontinuities, corresponding to spurious Fourier components at the larger wave vectors 
of the spectrum. Thereafter it is necessary to filter them out. A typical solution consists 
in multiplying the height profile h n by a windowing function that goes gradually to at 
the edge of the surface area and is almost constant around the centre; and then carrying 
on the Fast Fourier Transform with the adjusted data. This procedure is described in 
detail in Ref. [83]. 

Note also that from 



(h(x)h(o)) = J al 2 q C{q)e iq - X 

it follows that 



(h 2 ) = J d 2 q C(q). (D.2) 

But the mean of the square of the surface roughness amplitude can also be obtained 
directly from the experimental data via 

^ 2 ) = 4E(^-^) 2 - (D.3) 

n 

It is a good check of the numerical accuracy in the evaluation of the surface roughness 
power spectra to show that the (h 2 ) calculated from (D.2) and (D.3) gives identical 
results. 

Let us assume that the surface roughness is isotropic. In this case C(q) is 
independent of the direction of q, i.e., C(q) = C(q). We can make use of this fact 
to reduce the statistical noise in C(q) by performing an angular average in (D.2). Let 
us for simplicity write C(q) = Cirnqi) as C(m x ,m y ), where qi = 2n/L is the lattice 
spacing in the qr-space, m x = 0, 1, A — 1 and similarly for m y . In an analogous way 
we write C(q) = C(mq L ) as C(m), m = 1,2,..., A/2. We have 

Jd 2 qC(q) = 27rJ o 9l dqqC(q). (D.4) 
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The discretized version of (D.4) takes the form 

N/2 N/2 

9l E Yl C(m x ,m y ) = 2vrg| ^mC(m) (D.5) 

m x =0 m y =0 m 

where 

C(m x , m y ) = C(m x , m y ) + C(N - m x , m y ) 

+ C(m x , N - m y ) + C(N - m x , N - m y ) 

corresponds to C(q x , q y ) + C(—q x , q y ) + C(q x , —q y ) + C(—q x , —q y ) (remember that the 
Discrete Fourier Transform is periodic with period N). Eq. (D.5) is obeyed if we define 

CM = 2^1 E C(m> x ,m> y ) 

* n " l -l/2<\m'\~m<l/2 



With this definition the "sum-rule" 

N/2 

' '(III) 

m=l 



(h 2 ) = 2tt J dq qC{q) = 2nq 2 L ^ mC{i 
will be exactly obeyed. 




2 3 4 

logq (1/m) 



Figure Dl. Red line: calculated power spectrum for the surface shown in Fig. Al. 
Green line: the original power spectrum, employed to generate the surface. The surface 
extends over an area of 0.03 m x0.03m and the sampling length is a — 0.0001m. 



As an application in Fig. Dl we show the spectrum calculated with this method for 
the same surface described in Fig. Al. The agreement between the original spectrum 
and the calculated one provides evidence of the accuracy that can be achieved. 
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